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(54) Title: DYNAMICALLY VARIABLE FILTER 

(57) Abstract: A method of processing a signal pertaining to at least one electrical property of an organ of a subject is disclosed. 
The method comprises determining a physiological condition of the subject, selecting a frequency hand, filtering the signal acconling 
to the frequency band, and dynamically adapting the frequency band in response to a change in the physiological condition. 
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DYNAMICALLY VARL\BLE FILTER 

FIELD AND BACKGROUND OF THE INVENTION 

The present inventioa relates to processing of electrical signals, and more 
5 particularly to the filtering of a signal pertaining to at least one electrical property of an 
organ of a subject 

Technologies related to measurement of electrical i)roperties of organs, such as 
the measurement of bioimpedance are generally known. Typically, such technologies 
relate to the monitoring of physiological parameters by extracting physiologically 

10 significant characteristics from electrical measurements, see, e.g:, U.S. Patent No. 
6,577,897. Characteristics may include measures that aid in the discernment of 
physiological indications pertaining directly or indirectly to the state of organs (e.^., 
blood vessels, heart, lungs and the like), and reveal measures of various physiological 
conditions including critical life-threatening conditions. 

15 For example, heart diseases may be caused by (i) a failure in the autonomic nerve 

system v^ere the impulses from the central nervous system control to the heart muscle 
&il to provide a regular heart rate and/or (ii) an insufficient strength of the heart muscle 
itself where even though the patient has a regular heart rate, its force of contraction is 
insufficient. Either way, the amount of blood or the rate at uduch the blood is supplied 

20 by a diseased heart is abnormal and it is appreciated that an assessment of the state of a 
patient's circulation is of utmost importance. 

The simplest measurements, such as heart rate and blood pressure, may be 
adeqiiate for many patients, but if there is a cardiovascular abnormality then more 
detailed measurements are needed. 

25 Cardiac output (CO) is the volume of blood pumped by the heart during a time 

interval, which is typically taken to be a minute. Cardiac output is the product of heart 
rate (HR) and the amount of blood which is pumped with each heartbeat, also known as 
the stroke volume (SV). For example, the stroke volume at rest in the standing position 
averages between 60 and 80 ml of blood in most adults. Thus, at a resting heart rate of 

30 80 beats per minute the resting cardiac output varies between 4.8 and 6.4 L per min. 

A common clinical problem is that of hypotension (low blood pressure); this may 
occur because the cardiac output is low and/or because of low systemic vascular 
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resistance. This problem can occur in a wide range of patients, especially those in 
intensive care or postoperative high dependency units. In these high risk patients, more 
detailed monitoring is typically established including measuring central venous pressure 
via a central venous catheter and continuous display of arterial blood pressure via a 
5 peripheral arterial catheter. 

In addition to the above measurements, the measurement of cardiac output is 
useful. For example, when combined with arterial pressure measurements, cardiac 
ou^ut can be used for calculating die systemic vascular resistance. The measurement of 
cardiac output is useful both for establishing a padenf s initial cardiovascular state and 
10 for monitoring the response to various therapeutic interventions such as transfusion, 
infusion of inotropic drugs, infusion of vasoactive drugs (to increase or reduce systemic 
vascular resistance) or altering heart rate either pharmacologically or by adjusting pacing 
rate. 

Several methods of measuring cardiac output are presently known, representative 

15 Examples include the Pick method, described by Adolf Pick in 1870, the amount of 
oxygen taken up by the body during respiration and the difference in oxygen 
concentration between venous and arterial blood is used to calculate the cardiac output; 
the transoesophageal echocardiogmphy (see, e,g,^ U.S. Patent No. 6,142,941) in which 
cardiac output is derived £rom blood flow velocity (recorded via Doppler sMft) cross- 

20 sectional area of the blood vessel and heart rate; and the compliance based method (see, 
e.g,, U.S. Patent No. 6,485,431) in which the compliance of the arterial system is 
determined from measured arterial pressure and used for calculating the cardiac output 
as the product of the mean arterial pressure and compliance divided by a time constant 
Also known are catheter based methods such as thermodilution (see, e.g,^ U.S. Patent 

25 No. 4,153,048). 

A non-invasive method, known as thoracic electrical bioimpedance, was first 
disclosed in U.S. Patent No. 3,340,867 and has recently begun to attract medical and 
industrial attention [U.S. Patent Nos. 3,340,867, 4,450,527, 4,852,580, 4,870,578, 
4,953,556, 5,178,154, 5,309,917, 5,316,004, 5,505,209, 5,529,072, 5,503,157, 

30 5,469,859, 5,423,326, 5,685,316, 6,485,431, 6,496,732 and 6,511,438; 
U.S. Patent Application No. 20020193689]. The thoracic electrical bioimpedance 
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method has the advantages of providing continuous cardiac output measurement at no 
risk to the patient. 

Various methods employing bioimpedance are found in: International Patent 
Application Publication Nos. WO2004098376 and WO2006087696, U-S, Patent Nos. 
5 6,022,322, 5,615,689 and 5,642,734, and U.S. Published Application Nos. 
20030120170, 20060085048 and 20060122540, the contents of which are hereby 
incorporated by reference. 

SUMMARY OF THE INVENTION 
10 According to one aspect of embodunents the present invention there is provided 

a method of processing a signal pertaining to at least one electrical property of an organ 
of a subject. The method comprises determining a physiological condition of fiie 
subject, selecting a frequency band, filtering the signal according to the frequency band, 
and dynamically adapting the frequency band in response to a change in the 
15 physiological condition, thereby processing the signal. 

According to another aspect of embodiments of the present invention there is 
provided a filtering device. The filtering device comprises a first input unit for receiving 
an input pertaining to at least one electrical property of an organ of a subject, a second 
input unit for receiving data pertaining to a physiological condition of the subject, and a 
20 filtering unit configured for filtering the input signal accordmg to a frequency band 
which is dynamically ad^ted in response to a change in the physiological condition. 

According to yet another aspect of embodiments of the present invention there is 
provided a system for monitoring cardiac output, comprising the filtering device. 

According to still another aspect of embodiments of the present invention there is 
25 provided a system for predicting at least one of: a body cell mass, a fat free mass and 
total body water of a subject, comprising the filtering device. 

According to an additional aspect of embodiments of the present invention there 
is provided a system for determining hematocrit of blood in a body part of a subject, 
comprising the filtering device. 
30 According to yet an additional aspect of embodiments of the present invention 

there is provided a system for monitoring hydration status of a subject, comprising the 
filtering device. 
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According to still an additional aspect of embodiments of the present invention 
there is provided a system for discriminating tissue, comprising the filtering device. 

According to a further aspect of embodiments of the present invention there is 
provided a system for calculating the circumference of a body segment, comprising the 
filtering device. 

According to yet a further aspect of embodiments of the present invention there 
is provided a method of monitoring at least one electrical property of an organ of a 
subject The method comprises: sensing an input radiofirequency signal fix>m the organ, 
processing the input radiofirequency signal to provide a processed input signal, filtering 
the input signals using a dynamically variable filter to provide a filtered signal, and 
using the filtered signal for monitoring the at least one electrical property of the organ. 

According to still a further aspect of embodiments of the present invention there 
is provided apparatus for monitoring at least one electrical property of an organ of a 
subject. The apparatus comprises an input unit for receiving an input radiofirequency 
signal sensed firom the organ; a signal processing unit for processing the input 
radiofrequency signal to provide processed input signal; a filtering imit configured for 
filtering the input signal using dynamically variable filter to thereby provide a filtered 
signal; and a monitoring unit for monitoring the at least one electrical property of the 
organ based on the filtered signal. 

According to still a further aspect of embodiments of the present invention there 
is provided a system for monitoring at least one electrical property of an organ of a 
subject Hie system comprises a radiofirequency generator for generating an output 
radiofirequency signal and a plurality of electrodes, designed to be connectable to the 
skin of the subject, and configured for transmittmg the output radiofirequency signal to 
the organ and sensing an input radiofirequency signal firom the organ. The system fiuther 
comprises a monitoring apparatus, the apparatus described herein. 

According to further features in preferred embodiments of the invention 
described below, the dynamic2illy variable filter is adapted in response to a change in a 
physiological condition of the subject The filter is typically a band pass filter 
characterized by a firequency band defined by, e.g,^ a lower frequency bound and an 
upper firequency bound. 
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According to some embodiments of the present invention the physiological 
condition is a heart rate of the subject 

According to some embodiments of the present invention at least one of a lower 
bound of the frequency band and an upper bound of the frequency band parameter is a 
linear function of the heart mte. 

According to frirther features in preferred embodiments of the invention 
described below, a lower bound of the frequency band is about Fl(HR) where 
Fl(HR) = 0.9x(HR/60) Hz and HR is a heart rate of the subject in units of beats per 
minute. 

According to some embodiments of the present invention an upper boimd of the 
frequency band is about Fu(HR), \^diere Fu(HR) = 6 + L5x[(HR/60) - 1] Hz. 

According to some embodiments of tiie present invention the heart rate is 
determined from an ECG signal received from the subject. 

According to some embodiments of the present invention the upper frequency 
bound is determined using an iterative process. 

According to some embodiments of the present invention the iterative process is 
based on a comparison between a value of a physiological parameter as extracted from 
the filtered input signal and a value of the physiological parameter as extracted from a 
reference signal. According to some embodiments of the present invention the reference 
signal comprises the ECG signal. According to some embodiments of the present 
invention the physiological parameter is a ventricular ejection time (VET), 

According to some embodiments of the present invention each iteration of the 
iterative process comprises: if the comparison meets a predetermined criterion, then 
updating the upper frequency bound by calculating an average between a low threshold 
for the upper bound and a high threshold for the upper bound. The thresholds can be 
predetermined or they can be set in a previous iteration of the iterative process. 

According to some embodiments of the present invention the physiological 
According to some embodiments of the present invention the iterative process is 
terminated if a value of the VET as extracted from the filtered input signal is higher than 
a value of the VET as extracted from the reference signal. 

According to some embodiments of the present invention the VET is averaged 
over a plurality of heart beats. 
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According to some embodiments of the present invention tbs VET is extracted 
from an average heart beat morphology of the subject 

According to some embodiments of the present invention an initial value of the 
upper frequency bound m the iterative process is a linear fimction of the heart rate. 
According to some embodiments of tiie present invention the linear function is Fu(HR). 

According to some embodiments of the present invention the radiofrequency 
signal is filtered using an analog filter. 

According to some embodiments of the present mvention at least one quantity is 
calculated using the electrical property. The quantity can be a stroke volume, a cardiac 
output, a brain mtra luminal blood volume, a blood flow and the like. 

According to some embodiments of the present mvention the blood flow 
comprises at least one of: an external carotid blood flow rate, an internal carotid blood 
flow rate, an ulnar blood flow rate, a radial blood flow rate, a brachial blood flow rate, a 
common iliac blood flow rate, an extemal iliac blood flow rate, a posterior tibial blood 
flow rate, an anterior tibial blood flow rate, a peroneal blood flow rate, a lateral plantar 
blood flow rate, a medial plantar blood flow rate and a deep plantar blood flow rate. 

According to some embodiments of the present invention a phase shift of the 
input radiofrequency signal relative to an output radiofrequency signal transmitted to the 
organ is determined. The phase shift can be used for calculating the quantity or 
quantities. 

According to some embodiments of the present invention the amplitude 
modulation of the input radiofrequency signal is reduced or eUminated so as to provide 
signals of substantially constant envelope. 

According to some embodiments of the present invention a phase modulation of 
the input radiofrequency signal is maintained while reducing or eliminating the 
amplitude modulatioiL 

According to some embodiments of the present invention the input 
radiofrequency signal and the output radiofrequency signal are mixed so as to provide a 
mixed radiofrequency signal. According to some embodiments of the present invention 
the mixed radiofirequency signal comprises a radiofrequency sum and a radiofrequency 
difference. 
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According to some embodiments of the piesent invention the input 
radiofrequency signal is indicative of impedance the organ. According to some 
embodiments of the present invention the input radiofrequency signal is indicative of 
hemodynamic reactance of the organ. 

5 Embodiments of the present invention successfully address the shortcomings of 

the presently known configurations by providing techniques for processing a signal 
pertaining to one or more electrical property of an organ of a subject. 

Unless otherwise defined, all technical and scientific terms used herein have the 
same meaning as commonly understood by one of ordinary skill in the art to wMch this 

10 invention belongs. Although methods and materials similar or equivalent to those 
described herein can be xised in the practice or testii^ of the present invention, suitable 
methods and materials are described below. In case of conflict, the patent specification, 
including definitions, will control. In addition, the materials, methods, and examples are 
illustrative only and not intended to be limiting. 

15 Implementation of the method and system of the present invention involves 

performing or completing selected tasks or steps manually, automatically, or a 
combination thereof. Moreover, according to actual instrumentation and equipment of 
preferred embodiments of the method and system of the present invention, several 
selected steps could be implemented by hardware or by software on any operating 

20 system of any firmware or a combination thereof. For example, as hardware, selected 
steps of the invention could be implemented as a chip or a circuit As software, selected 
steps of the invention could be implemented as a plurality of software instructions being 
executed by a computer using any suitable operating system. In any case, selected steps 
of the method and system of the invention could be described as being performed by a 

25 data processor, such as a computing platform for executing a plurality of instructions. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the 
accompanying drawings. With specific reference now to the drawings in detail, it is 
30 stressed that the particulars shown are by way of example and for purposes of illustrative 
discussion of the preferred embodiments of the present invention only, and are presented 
in the cause of providing what is believed to be the most useful and readily imderstood 
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description of the principles and conceptual aspects of the invention. In this regard, no 
attempt is made to show structural details of the invention in more detail than is 
necessary for a fundamental understanding of the invention, the description taken with 
the drawings making apparent to those skilled in the art how the several forms of the 
5 invention may be embodied in practice. 
In the drawings: 

FIG. 1 is a schematic illustration of a filtering device, according to various 
exemplary embodiments of the present invention; 

FIGs, 2a-b show a representative example of dynamically varying frequency 
10 bounds, employed according to embodiments of the present invention; 

FIG. 2c show a representative example of a dynamically varying frequency band, 
employed according to embodiments of the present invention; 

FIG. 3a is a flowchart diagram of a method suitable for processing an input 
signal, according to various exemplary embodiments of the present invention; 
15 FIG. 3b is a flowchart diagram of an iterative process for selecting a frequency 

bound, according to various exemplary embodiments of the present invention; 

FIG. 3c is a schematic illustration of a procedure for extracting ventricular 
ejection time from an EGG signal; 

FIG. 3d is a schematic illustration of a procedure for extracting ventricular 
20 ejection time from a derivative of a filtered signal according to various exemplary 
embodiments of the present invoition; 

FIGs. 3e-f are schematic illustrations of beat morphology characteri2:ation 
procediures, according to various exemplary embodiments of the present invention; 

FIG. 4 is a schematic illustration of apparatus for monitoring one or more 
25 electrical properties of an organ of a subject, according to various exemplary 
embodiments of the present invention; 

FIG. 5 is a schematic illustration of a signal processing unit, according to various 
exemplary embodiments of the present invention; 

FIG, 6 is a block diagram of electronic circuitry, according to various exemplary 
30 embodiments of the present invention; 
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FIG. 7 is a schematic illustration of a system for monitoring at least one electrical 
property of an organ of a subject, according to a preferred embodiment of the present 
invention; 

FIGs. 8a-e are schematic illustrations showing perspective (Figure Sa), ftont 
5 (Figure 8bX rear (Figure 8c), side (Figure 8d) and top (Figure 8e) views of a sticker 
which can be used for transmitting and sensing the ladiofrequency signal, according to 
one embodiment of the present invention; 

FIG. 8f is a schematic illustration of a package of several stickers accorcUng to an 
embodiment of the present invention; 
10 FIGs. 9a-d are schematic illustrations of various electronic circuitries, according 

to exemplary embodiments of the present invention; 

nGs. lOa-lO-e, lla-lle, 12a-12e, 13a-13e, 14a.l4e, 15a-15e, 16a-16g and 17a- 
17g show snapshots of the display of a prototype system, manufactured and configured 
according to various exemplary embodiments of the present invention; and 
15 FIGs. 18a-18b, 19a-19b and 20a-20b are plots of cardiac output as calculated 

from a signal filtered according to various exemplary embodiments of the present 
invention. 

DESCRIPTION OF EXEMPLARY EMBODIMENTS 

20 The present embodiments comprise a method, device apparatus and system 

which can be used for processing a signal. Specifically, but not exclxisively, the present 
embodiments can be used for processing a radiofrequency signal sensed from an organ 
of a subject and for monitoring one or more electrical properties of an organ, e,g, , for the 
purpose of determining one or more quantities which are related to electrical properties. 

25 Tlius, for example, exemplary embodiments of the present invention can be used for 
calculating stroke volume, cardiac output, brain intra luminal blood volume and/or blood 
flow. Embodiments of the present invention can also be used for discriminating tissue 
and/or determining at least one of: body cell mass, fat free mass, total body water, 
hematocrit of blood, hydration status and circumference of a body segment. 

30 The principles and operation of a method, device apparatus and system according 

to the present embodiments may be better understood witti reference to the drawings and 
accompanying descriptions. 
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Before explaining at least one embodiment of the invention in detail, it is to be 
understood that the invention is not limited in its application to the details of 
construction and the arrangement of the components set fortii in the folloAving 
description or illustrated in the drawings. The invention is capable of oth^ 
5 embodiments or of being practiced or carried out in various ways. Also, it is to be 
imderstood that the phraseology and terminology employed herein is for the purpose of 
description and should not be regarded as limiting. 

Computer programs implementing the method according to embodiments of the 
present invention can commonly be distributed to users on a distribution medium such 

10 as, but not limited to, a floppy disk, CD-ROM and flash memory cards. From the 
distribution medium, the computer programs can be copied to a hard disk or a similar 
intermediate storage medium. The computer programs can be run by loading the 
computer instructions either from their distribution medium or their intermediate storage 
medium into the execution memory of the computer, configuring the computer to act in 

15 accordance with llie method of this invention. All these operations are well-known to 
those skilled in the art of computer systems. 

A typical system for monitoring electrical properties of a body section, such as a 
bioimpedance system, includes a tetrapolar array of circumferential band electrodes 
cormected to the subject at the base of the neck and surrounding the circumference of the 

20 lower chest, at the level of the xiphoid process. When a constant magnitude alternating 
current flows through the upper cervical and lower thoracic band electrodes, a voltage, 
proportional to the thoracic electrical impedance (or reciprocally proportional to the 
admittance), is measured between the inner cervical and thoracic band electrodes. The 
portion of the cardiac synchronous impedance change, temporally concordant with the 

25 stroke volume, is ascribed solely and uniquely to volume changes of the aorta during 
expansion and contraction over the heart cycle. A typical printed circuit board of such 
system comprises one or more band pass filters, a half-wave rectification circuit and one 
or more low pass filters. 

The present Inventor discovered techniques for reducing the noise associated 

30 with conventional systems. As demonstrated in the Examples section that follows, the 
present Inventor succeeded in reducing noise introduced due to patient agitation or other 
physiological phenomena like breathing. The present Inventor discovered techniques for 
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separating and differentiating between a cardiovascular bioreactance signal and a 
respiratory bioreactance signal, where the latter is typically much larger than the former. 

The present Inventor has realized that the noise level is proportional to the 
bandwidth of the band pass filter and that a considerable portion of the noise passes the 

5 band pass filter hence being folded into the half-wave rectification circuit 

The present Inventor also discovered techniques for reducing or eliminating AM 
noise hence significantly improving the ability to provide accurate measurement 

Referring now to the drawings. Figure 1 illustrates a filtering device 10, 
according to various exemplary embodiments of the present invention. Device 10 

10 comprises a first input unit 12 which receives an input signal 16 pertaining to one or 
more electrical properties of an organ of a subject For example, signal 16 can relate to 
tibie hemodynamic reactance of the organ. 

As used herein, "hemodynamic reactance" refers to the imaginary part of the 
impedance. Techniques for extracting the imaginary part fix>m the total impedance are 

15 known in the art. Typically, such extraction is performed at hardware level but the use 
of algorithm at a software level is not excluded from the scope of the present invention. 
Signal 16 can be provided, for example, by processing a radiofi'equency signal sensed 
fi:om the organ, as fiirther detailed hereinunder. 

In various exemplary embodiments of the invention device 10 fiirther comprises 

20 a second input unit 14 which receives data 18 pertaining to a physiological condition of 
the subject The physiological condition is preferably, but not obligatorily, the heart rate 
of the subject, and the data pertaining to the physiological condition can be analog data 
or digital data, as desired. While the embodiments below are described with a particular 
emphasis to physiological condition which is a heart rate, it is to be understood that 

25 more detailed reference to the heart rate is not to be interpreted as limiting the scope of 
the invention in any way. For example, in exemplary embodiments of the present 
invention the physiological condition is a ventilation rate of the subject, a repetition rate 
of a particular muscle unit and/or one or more characteristics of an action potential 
sensed electromyography. 

30 Device 10 further comprises a filtering unit 20 which filters the input signal 16 to 

provide a filtered signal 22. In various exemplary embodiments of the invention the 
filtering is according to a frequency band which is dynamically adapted in response to a 
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change in the physiological condition of the subject. It was found by the Inventor of the 
present invention that the dynamical adaptation of the frequency band to tiie 
physiological condition of the subject can significantly reduce the influence of unrelated 
signals on the measured or monitorii^ of electrical properties of the body section. 
5 The adaptation of the frequency band to the physiological condition can be 

according to any adaptation scheme known in the art. For example, one or more 
parameters of the frequency band (e.g:, lower bound, upper bound, bandwidth, central 
frequency) can be a linear function of a parameter characterizing the physiological 
condition. Such parameter can be, for example, the number of heart beats per minute. 

10 A representative example of a dynamically varying frequency bounds, employed 

according to some embodiments of the present invention by unit 20, is illustrated in 
Figures 2a-b. Shown in Figures 2a-b is the functional dependence of the frequency 
bounds (upper bound in Figure 2a and lower bound in Figure 2b) on the heart rate of the 
subject. As shown in Figure 2a, the upper bound of the frequency band varies linearly 

15 such that at a heart rate of about 60 beats per minute (bpm) the upper bound is about 6 
Hz, and at a heart rate of about 180 bpm the upper bound is about 9 Hz. As shown in 
Figure 2b, the lower bound of the frequency band varies linearly such that at a heart rate 
of about 60 the lower boimd is about 0.9 Hz bpm and at a heart rate of about 180 bpm 
the lower bound is about 2,7 Hz. 

20 As used herein the term "about" or "approximately" refers to ± 10 %. 

In some embodiments of the present invention the upper bound approximately 
equals the function Fu(HR) defined as Fu(HR) = 6 + 1.5 x [(HR/60) - l]Hz, where HR is 
the heart rate of the subject in units of bpm. In some embodiments, the upper bound 
equals Fu(HR) at all times, while in other embodiments, the upper bound is set using an 

25 iterative process. 

In some embodiments of the present invention the lower bound approximately 
equals the function Fl(HR) defined as Fl(HR) = 0.9 x (HR/60) Hz. In some 
embodiments, the lower bound equals Fl(HR) at all times while in other embodiments 
the lower bound is set by an iterative process. 

30 Representative examples of iterative process suitable for some embodiments of 

the present invention are provided hereinunder. 
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A dynamically varsdng band pass filter (BPF) characterized by a dynamically 
var^g upper frequency bound and a dynamically varying lower firequency bound, 
according to some embodiments of the present invention is illustrated in Figure 2c. As 
shown, each heart rate is associated with a frequency band defined by a lower bound and 

5 an upper bound. For example, for a heart rate of 60 bpm. Figure 2c depicts a BPF in 
which the lower bound is about 0.9 Hz and the upper bound is about 6 Hz. 

It is to be understood that the values presented above and the functional relations 
illustrated in Figures 2a-b are exemplary embodiments and should not be considered as 
limiting the scope of the present invention in any way. In other exemplary 

10 embodiments, the fimctional relations between the frequency band and the physiological 
condition can have different slopes and/or ofisets, or they can be non-linear. 

Figure 3a is a flowchart diagram of a method suitable for processing an input 
signal, according to various exemplary embodiments of the present invention. The 
method can be executed by activating filtering device^ 10 or using any other filtering 

15 device supplemented by appropriate circuitry. Selected steps of the method can be 
embodied in many forms. For example, the selected steps can be embodied in on a 
tangible medium such as a computer for performing the selected steps. The selected 
steps can be embodied on a computer readable medium, comprising computer readable 
instructions for carrying out the selected steps. The selected steps can also be embodied 

20 in electronic device having digital computer capabilities arranged to run the computer 
program on the tangible medium or execute the instruction on a computer readable 
medium. 

The method begins at step 30 and continues to step 31 in which the physiological 
condition of the subject is determined. The physiological condition can be, as stated, a 

25 heart rate and it can be determined using any procedure known in the art, such as, but 
not limited to, analysis of an electrocardiogram (ECG) signal or the like. The method 
continues to step 32 in which a frequency band is selected based on the physiological 
condition of the subject, and proceeds to step 33 in N\iiich the input signal is filtered 
according to frequency band. In various exemplary embodiments of the invention the 

30 method loops back to step 31 so as to dynamically adapt the frequency band in response 
to a change in the physiological condition. 
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The selection of frequency band can be according to any adaptation scheme, 
including, without limitation, the use of one or more linear functions (e.^., the functions 
Fu and Fl) as further detailed hereinabove, and/or an it^tive process as further detailed 
hereinbelow. 
5 The method ends at step 34. 

Following is a description of an iterative process for determining the frequency 
band of the band pass filter which filters the input signal according to some 
embodiments of the present invention. The iterative process can, in some embodiments, 
based a comparison between a value of a physiological parameter as extracted or 
10 calculated &om the filtered input signal and a value of the same physiological parameter 
as extmcted or calculated from a reference signal, for example, an BCG signal. 

Hie term physiological parameter refers to any variable parameter which is 
measurable or calculable and is representative of a physiological activity, particularly, 
but not necessarily, activity of tiie heart In various exemplary embodiments of the 
15 invention ttie physiological parameter is other than the heart rate per se. The 
physiological parameter can be a time-related parameter, amplitude-related parameters 
or combination thereof. 

Typically, the filter sign^ and the reference signal are expressed in terms of 
amplitude as a function of the time. Thus, time-related parameters are typically 
20 calculated using abscissa values of the signals and amplitude-related parameters are is 
typically calculated using ordinate values of the signals. 

Representative of time-related physiological parameters suitable for the present 
embodiments include, without limitation, systolic time, diastolic time, pre-ejection 
period and ejection time. A representative example of amplitude-related physiological 
25 parameter suitable for the present embodiments includes, without limitation, cardiac 
contractility, maximal amplitude above zero during a single beat, maximal peak-to-peak 
amplitude during a single beat, and RMS level during a single beat Also contemplated 
are various slopes parameters, such as, but not limited to, the average slope between two 
points over the signal. 

30 In various exemplary embodiments of the invention the physiological parameter 

is a ventricular ejection time (VET). 
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While the embodiments below are described with a particular emphasis to VET 
as the physiological parameter, it is to be understood that more detailed reference to 
VET is not to be inteipreted as limiting the scope of the invention in any way. 

The present inventors discovered that a significant amount of the biological 
S information for a particular subject can be obtmned firom a frequency range between 
Fl(HR) and S.S Hz, where HR is the heart rate of the subject. It was further discovered 
by the present inventors that for some medical conditions some of the information can 
reside between 5.5 Hz and Fu(HR). 

The advantage of the comparison between two different techniques for extracting 
10 or calculating the same physiological parameter, is that it allows to substantially 
optimize the upper firequency bound of the of the band pass filter. In various exemplaiy 
embodiments of the invention in each iteration of the iterative process,, the comparison is 
repeated. If the comparison meets a predetermined criterion, the upper frequency bound 
is updated by calculating an average between a low threshold for the upper bound and a 
15 high threshold for the upper bound. The lower frequency boimd can be a constant 
bound, e,g.^ a constant frequency ^^ch is from about 0.9 Hz to about 2.7 Hz), or it can 
be dynamic, e.g., Fl(HR), HR being the heart rate of the subject before or during the 
respective iteration. 

Hie low and high thresholds for the upper bound can be set in more than one 
20 way. In some embodiments, the low and high thresholds are predet^mined (namely 
they determined a priori before the iterative process), in some embodiments, the 
thresholds are set in a previous iteration of iterative process, in some embodiments one 
of the thresholds is predetermined and the other threshold is set in a previous iteration of 
iterative process. In any event, the first iteration is based on two thresholds which are 
25 determined a priori before the iterative process. It was found by the inventors of the 
present invention that, at least initially (/.e., at the first iteration), the first threshold can 
be about Fu(40), which in various exemplary embodiments of the invention is about 
5.5 Hz, and the second threshold can be the calculated value of Fu(HR), HR being the 
heart rate of the subject before or during the respective iteration. 
30 The predetermined criterion used during the iterations can be, for example, that 

the results of the two calculations are similar {e.g, within about 40 % or 30 % or 25 % 
of each other). The predetemiined criterion can also relate to the direction of difference 
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between the two calculations. Broadly, for time-related parameters, the upper bound is 
updated if the value of the parameter as calculated based on the reference signal is 
higher than value of the parameter as calculated based on the filtered signal, and for 
amplitude-related parameters the upper bound is updated if the value of parameter as 
5 calculated based on the reference signal is lower than the value of the parameter as 
calculated based on the filtered signal. For slope-related parameters, the upper bound is 
typically updated if the value of the parameter as calculated based on the reference 
signal is higher than the value of the parameter as calculated based on the filtered signal. 
A Boolean combination between the above criteria can also be used as a 

10 criterion. For example, an AND Boolean combination can be enq)loyed in which case 
the upper firequency bound can be updated if the results of the two calculations are 
similar and the calculation according to the filtered signal indicates an abnormal 
physiological condition while the calculation according to the reference signal indicates 
a normal physiological condition. 

15 Figure 3b is a flowchart diagram of an iterative process for selecting the upper 

firequency bound, according to various exemplary embodiments of the present invention. 
The description is for a physiological parameter which is VET, but, as stated, it is not 
intended to limit the scope of the present invention to this type of physiological 
parameter. 

20 The iterative process begins at 60 and continues to 61 in which the upper 

firequency bound is set to the value of Fu(HR), HR being the heart rate of the subject 
before the initiation of the iterative process. The heart rate can be inputted or it can be 
determined by the process, e.g,, from the ECG signal. The process continues to 62 in 
which initial values are assigned to two frequency thresholds. The frequency thresholds 

25 are denoted in Figure 3b by Tl and T2. In various exemplary embodiments of the 
invention the initial value of Tl is Fu(40) and the initial value of T2 is the initial upper 
frequency bound. 

The process continues to 63 and 64 in which VET is calculated separately from 
the filtered input signal (63) and from a reference signal (64), e.g., ECG. The VET as 
30 calculated from the input signal is denoted in Figure 3b VETl and the VET as calculated 
from the ECG signal is denoted in Figure 3b VET2. In the first iteration, the filtered 
signal from which VETl is calculated is preferably obtained by filtering the input signal 
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using a band pass filter defibaed between the upper and lower frequency bounds as 
initially set at 61. 

The process continues to decision 65 in vMch VETl and VET2 are compared. If 
VETl is higher than VET2, the iterative process continues to 73 where it is terminated. 
5 If VETl is not higher than VET2, the process continues to 66 in \Auch the upper boimd 
is updated to the average AVE of Tl and T2. In various exemplary embodiments of the 
mvention AVE(T1, T2) is the arithmetic mean of Tl and T2 AVE(T1, 
T2) = (Tl+T2)/2), this need not necessarily be the case, since, in some embodiments, a 
different averaging scheme (cg.-^ a weighted average, a geometric mean, harmonic 
10 mean, RMS, etc.) can be employed. 

From 66 the process continues to 67 in whidi VETl is recalculated but from a 
signal which is filtered using an updated band pass filter. The lower bound of the 
iqxlated band pass filter can be the initial lower boimd or it can be updated based on the 
heart rate of the subject immediately before tfie filtration of the input signal (e.g,, 
15 according to Fl described above). The upper bound of the updated band pass filter is 
preferably the updated upper bound. Optionally, Ihe process continues to 68 m which 
the VET2 is also recalculated ftom the reference signal. Alternatively, the value of 
VET2 fi:om 64 can be used. 

The process then continues to decision 69 in which the process, determines 
20 whether or not a predetermined termination criterion is met If the predetermined 
termination oiterion is met, the iterative process continues to 73 where it is terminated, 
otherwise the process continues to decision 70 in which VETl and VET2 are compared. 
If the deviation between VETl and VET2 is lower than a predetermined threshold A and 
VETl is lower than VET2, the process continues to 71 at v^ch the value of the upper 
25 bound is assigned to T2, otherwise the process contmues to 72 at which the value of the 
upper bound is assigned to Tl. 

From 71 and 72 the process loops back to 66 at which an additional iteration 

begins. 

The threshold A employed at decision 70 is typically expressed as a fraction of 
30 VET2 or a fraction of (VET2.VET1)A/^ET2. In various exemplary embodiments of the 
invention A = p* VET2, where p < 0.5, for example, p 0.4 or p = 0.3 or p = 0.25. 
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The termination criterion employed at decision 69 can be for example, a maximal 
number of iterations. In this embodiment, the process counts the number of iterations 
and compares them to a predetermined iteration number threshold. If the number of 
iterations exceeds iteration ntmiber threshold the process detennines that the termination 
5 condition is met and continues to end 73. Typically, but not necessarily, the value of the 
predetermined iteration number threshold is from about 3 iterations to about 10 
iterations, e.g., 5 or 6 iterations. 

The termination criterion can also include a comparison of the value of VETl to 
a predetermined VET threshold Tvet, which may be absolute, subject-specific or relative 

10 to VET2. It this embodiment, the process compares the value of VETl to Tvet- If 
VETl is higher than Tvet the process determines that the termination condition is met 
and continues to end 73. Typically, but not necessarily, a relative value of Tvet is 
employed. For example, in various exemplary embodiments of the invention 
Tvet = p*VET2 where p < 0.5, for example, p = 0.4 or p = 0.3 or p = 0.25. 

15 In any of the above embodiment, the VET (either VETl or VET2) can be 

calculated over a single heart beat, or, more preferably, it can be averaged over two or 
more heart beats. The advantage of using an averaging procedure rather than using a 
single beat is that it attenuates random disturbances wiiich may be present in the signal 
during the iterative process. Nevertheless, extraction of the VET from a single beat is 

20 not excluded from the scope of the present invention. 

For a single heart beat, the calculation of VET can be performed by 
characterizing the morphology of the beat, identifying two or more identifiable points on 
the beat and measuring the time between the identified points. 

Figure 3c illustrates a procedure for extracting VET2 when the reference signal 

25 is an ECG signal. Shown in Figure 3c is a typical morphology of a single beat of an 
ECG signal as a function of the time. VET2 can be defined as the time period 
(difference between the abscissa values) between the R peak and the T peak of the ECG 
signal. 

When the filtered signal is hemodynamic reactance, the value of VETl is 
30 preferably extracted from the first derivative of the filtered signal. The procedure is 
illustrated in Figure 3d, which illustrates a typical morphology of a single beat of the 
hemodynamic reactance N and its first derivative dW, as a function of the time. As 



wo 2009/022330 



PCT/IL2008/001105 



19 

shown, dN has two zeroes Oi and O2 over the beat, with a pomt of local maximum Mi 
between the zeroes and a point of local minimum M2 after the second zero. In some 
embodiments of the present invention VETl is defined as the time period (difference 
between the abscissa values) between the first zero Oi and the first minjymim M2 after 
5 the second zero O2. 

The beat motphology of the filtered signal can be characterized by identifying 
identifiable points on the reference signal and using the tune associated with these pomts 
(abscissa values) for defining anchor points on the filtered signal. 

Two exemplary beat morphology characterization procedure of the filtered signal 
10 for embodiments in which the input signal is the hemodynamic reactance and the 
reference signal is the ECG signal are illustrated in Figures 3e-f. 

In Figure 3e, a single beat of dN (first derivative of the hemodynamic reactance 
N) is defined between two anchor endpoints: a first (left) endpoint has the abscissa value 
of the Q peak of the ECG signal, and a second (right) endpoint has the abscissa value of 
15 the R peak of the next ECG signal. In other words, a single beat of diV^ has a width 
wiiich equals the following sum of five successive intervals over the ECG: 
QR+RS + ST + TQ + QR. 

In Figure 3f, a single beat of dATis defined using three anchor points: the two 
anchor endpoints as described in Figure Be and an intermediate anchor point which has 
20 the abscissa value of the global maximum A of the hemodynamic reactance N between 
the two endpoints. 

When the VET is averaged, the calculation can be done in more than one way. 
In some embodiments, the same morphology characterization is employed for a 
plurality of beats over a predetermined interval of the respective signal, so as to provide 
25 one local VET for each beat. The VET can be defined as the average of all local VETs. 
Any averaging procedure can be employed, include, without limitation, arithmetic mean, 
weighted average, geometric mean, harmonic mean, RMS and the like. 

In some embodiments, the morphologies of the beats are averaged over an 
ensemble of beats in the filtered signal, to provide an average beat morphology, and the 
30 VET is be determined by identifying two or more identifiable points on the average beat 
morphology and measuring the time between the identified points. For example, for 
each beat in the ensemble, the morphology can be characterized as described above and 
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all tbie morphologies can be averaged, e.g., point by point. The morpholo^es can also 
be averaged in segments. This embodiments is particularly useful when a single beat is 
defined using more than two anchor points in which case each segment over the beats 
(between two successive anchor points) can be aveiaged, e,g., point-by-point The 
5 average beat morphology can then be obtained by stitching the avemged segments. For 
example, in the embodiment illustrated in Figure 3f the beat is defined using two 
intermediate anchor points and an intermediate point In tibds embodiment, each beat has 
a left segment (fix)m the abscissa value of Q to the abscissa value of A) and a right 
segment (fix>m the abscissa value of A to the abscissa value of R). The left segments of 

10 all beats in the ensemble can be averaged to provide a left segment average, and the right 
segments of all beats in the ensemble can be averaged to provide a right segment 
average. The average beat morphology can be obtained by stitching the left segment 
average to the right segment average. 

In various exemplary embodiments of the invention the time scale of the beats in 

15 the ensemble is adjusted so as to fit all the beats in the ensemble to a single time scale. 

The result of this average is a single beat morphology fi:om which VETl can be 
extracted as described above. For example, when the first derivative of the 
hemodynamic reactance is used, the average morphology typically has the shape 
illustrated in Figure 3d. VETl can then be extracted as the time period between Oi and 

20 02- The first derivative of the signal can be calculated before or after averaging. When 
the first derivative is calculated before the averaging, the averaging procedure described 
above is performed with respect to the derivative of the signal. When the fu^ derivative 
is calculated after the averaging, the averaging procedure described above is performed 
with respect to the signal and the obtained average is then differentiated. Calculation of 

25 the first derivative after averaging is preferred from the standpoint of noise reduction. 

The predetermined time period over which an averaging procedure is performed 
to extract the VET typically extends over about 10 heart beats. 

A particular advantage of the device and method of the present embodiments is 
that they can be implemented in many systems designed for measuring or monitoring 

30 electrical properties of body sections, thereby improving their performance, e.g.^ by 
increasing their signal to noise ratio at least for situations in which the amount of noise 
is high. Representative examples of such systems include, without limitation, a system 
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DYNAMICALLY VARIABLE FILTER 

FIELD AND BACKGROUND OF THE INVENTION 

The present invention relates to processing of electrical signals, and more 
5 particularly to the filtering of a signal pertaining to at least one eiectiical property of an 
organ of a subject 

Technologies related to measurement of electrical properties of organs, such as 
the measurement of bioimpedance are generally known. Typically, such technologies 
relate to the monitoring of physiological parameters by extracting physiologically 

10 significant characteristics from electrical measurements, see, e.g:, U.S. Patent No. 
6,577,897. Characteristics may include measures that aid in the discernment of 
physiological indications pertaining directly or indirectly to the state of organs (e.g., 
blood vessels, heart, lungs and the like), and reveal measures of various physiological 
conditions including critical life-threatening conditions. 

15 For example, heart diseases may be caused by (i) a failure in the autonomic nerve 

system where the impulses firom the central nervous system control to the heart muscle 
&il to provide a regular heart rate and/or (ii) an insufficient strength of the heart muscle 
itself where even though the patient has a regular heart rate, its force of contraction is 
insufficient Either way, the amount of blood or the rate at which the blood is supplied 

20 by a diseased heart is abnormal and it is appreciated that an assessment of the state of a 
patient's circulation is of utmost importance. 

The simplest measurements, such as heart rate and blood pressure, may be 
adequate for many patients, but if there is a cardiovascular abnormality then more 
detailed measurements are needed. 

25 Cardiac output (CO) is the volume of blood pumped by the heart during a time 

interval, which is typically taken to be a minute. Cardiac output is the product of heart 
rate (HR) and the amount of blood which is pumped vnttx each heartbeat, also known as 
flie stroke voliune (SV). For example, the stroke volume at rest in the standing position 
averages between 60 and 80 ml of blood in most adults. Thus, at a resting heart rate of 

30 80 beats per minute the resting cardiac output varies between 4.8 and 6.4 L per min. 

A common clinical problem is that of hypotension (low blood pressure); this may 
occur because the cardiac output is low and/or because of low systemic vascular 
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DYNAMICALLY VARIABLE FILTER 

FIELD AND BACKGROUND OF THE INVENTIQN 

The present invention relates to processing of electrical signals, and more 
5 particularly to the filtering of a signal pertaining to at least one electrical property of an 
organ of a subject 

Technologies related to measurement of electrical properties of organs, such as 
the measurement of bioimpedance are generally known. Typically, such technologies 
relate to the monitoring of physiological parameters by extracting physiologically 

10 significant characteristics from electrical measurements, see, e.g.^ U.S. Patent No. 
6,577,897. Characteristics may include measures that aid in the discernment of 
physiological indications pertaining directiy or indirectly to the state of organs (e.^., 
blood vessels, heart, lungs and the like), and reveal measures of various physiological 
conditions including critical life-threatening conditions. 

15 For example, heart diseases may be caused by (i) a failure in the autonomic nerve 

system where the impulses from the central nervous system control to the heart muscle 
faH to provide a regular heart rate and/or (ii) an insufficient strength of the heart muscle 
itself where even though the patient has a regular heart rate, its force of contraction is 
insufficient Either way, the amount of blood or the rate at Mdiich the blood is supplied 

20 by a diseased heart is abnormal and it is appreciated that an assessment of tiie state of a 
patienfs circulation is of utmost importance. 

The simplest measurements, such as heart rate and blood pressure, may be 
adequate for many patients, but if there is a cardiovascular abnormality then more 
detailed measurements are needed. 

25 Cardiac output (CO) is the volume of blood pumped by the heart during a time 

interval, which is typically taken to be a minute. Cardiac output is the product of heart 
rate (HR) and the amount of blood which is pumped with each heartbeat, also known as 
the stroke volume (SV). For example, the stroke volume at rest in the standing position 
averages between 60 and 80 ml of blood in most adults. Thus, at a restmg heart rate of 

30 80 beats per minute the resting cardiac output varies between 4.8 and 6.4 L per min. 

A common clinical problem is that of hypotension (low blood pressure); this may 
occur because the cardiac output is low and/or because of low systemic vascular 
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resistance. This problem can occur in a wide range of patients, especially those in 
intensive care or postoperative high dependency units. In these high risk patients, more 
detailed monitoring is typically established including measuring central venous pressure 
via a central venous catheter and continuous display of arterial blood pressure via a 
5 peripheral arterial catheter. 

In addition to the above measurements, the measurement of cardiac output is 
usefiil. For example, when combined with arterial pressure measurements, cardiac 
output can be used for calculating die systemic vascular resistance. The measurement of 
cardiac output is useM both for establishing a patient's initial cardiovascular state and 
10 for monitoring the response to various therapeutic interventions such as transfusion, 
infusion of inotropic dmgs, infusion of vasoactive drugs (to increase or reduce systemic 
vascular resistance) or altering heart rate either pharmacologically or by adjusting pacing 
rate. 

Several methods of measuring cardiac output are presently known, representative 

15 Examples include the Pick method, described by Adolf Pick in 1870, the amount of 
oxygen taken up by the body during respiration and the difference in oxygen 
concentration between venous and arterial blood is used to calculate the cardiac output; 
the transoesophageal echocardiography (see, e.g.^ U.S. Patent No. 6,142,941) in which 
cardiac output is derived from blood flow velocity (recorded via Doppler shift) cross- 

20 sectional area of the blood vessel and heart rate; and the compliance based method (see, 
e.g., U.S. Patent No. 6,485,431) in which the compliance of the arterial system is 
detemiined from measured arterial pressure and used for calculating the cardiac output 
as the product of the mean arterial pressure and compliance divided by a time constant 
Also known are catheter based methods such as thermodilution (see, e.g., U.S. Patent 

25 No. 4,153,048). 

A non-invasive method, known as thoracic electrical bioimpedance, was first 
disclosed in U.S. Patent No. 3,340,867 and has recently begun to attract medical and 
industrial attention [U.S. Patent Nos. 3,340,867, 4,450,527, 4,852,580, 4,870,578, 
4,953,556, 5,178,154, 5,309,917, 5,316,004, 5,505,209, 5,529,072, 5,503,157, 

30 5,469,859, 5,423,326, 5,685,316, 6,485,431, 6,496,732 and 6,511,438; 
U.S. Patent Application No. 20020193689]. The thoracic electrical bioimpedance 
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method has the advantages of providing continuous cardiac output measurement at no 
risk to the patient 

Various methods employing bioimpedance are found in: International Patent 
Application Publication Nos. WO2004098376 and WO2006087696, U-S. Patent Nos. 
5 6,022,322, 5,615,689 and 5,642,734, and U.S. Published Application Nos. 
20030120170, 20060085048 and 20060122540, the contents of which are hereby 
incorporated by reference. 

SUMMARY OF THE INVENTION 

10 According to one aspect of embodiments the present invention there is provided 

a method of processing a signal pertaining to at least one electrical property of an organ 
of a subject The method comprises determining a physiological condition of the 
subject, selecting a frequency band, filtering the signal according to the frequency band, 
and dynamically ad^ting the frequency band in response to a change in ttxe 

15 physiological condition, thereby processing the signal. 

According to another aspect of embodiments of the present invention there is 
provided a filtering device. Hie filtering device comprises a first input unit for receiving 
an input pertsdning to at least one electrical property of an organ of a subject, a second 
input unit for receiving data pertaining to a physiological condition of the subject, and a 

20 filtering unit configured for filtering the input signal according to a frequency band 
which is dynamically adapted in response to a change in the physiological condition. 

According to yet another aspect of embodiments of the present invention there is 
provided a system for monitoring cardiac output, comprising the filtering device. 

According to still another aspect of embodiments of the present invention there is 

25 provided a system for predicting at least one of: a body cell mass, a fat free mass and 
total body water of a subject, comprising the filtering device. 

According to an additional aspect of embodiments of the present invention there 
is provided a system for determining hematocrit of blood in a body part of a subject, 
comprising the filtering device, 

30 According to yet an additional aspect of embodiments of the present invention 

there is provided a system for monitoring hydration status of a subject, comprising the 
filtering device. 
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According to still an additional aspect of embodiments of the present invention 
there is provided a system for discriminating tissue, comprising the filtering device. 

According to a further aspect of embodiments of the present invention there is 
provided a system for calculating the circumference of a body segment, comprising the 
filtering device. 

According to yet a further aspect of embodiments of the present invention there 
is provided a method of monitoring at least one electrical property of an organ of a 
subject. The method comprises: sensing an input radiofrequmcy signal from the organ, 
processing the input radiofirequency signal to provide a processed input signal, filtering 
the input signals using a dynamically variable filter to provide a filtered signal, and 
usmg the filtered signal for monitoring the at least one electrical property of the organ. 

According to still a further aspect of embodiments of the present invention there 
is provided apparatus for monitoring at least one electrical property of an organ of a 
subject. The apparatus comprises an input unit for receiving an input radiofi:equency 
signal sensed from the organ; a signal processing unit for processing the input 
radiofrequency signal to provide processed input signal; a filtering unit configured for 
filtermg the input signal using dynamically variable filter to thereby provide a filtered 
signal; and a monitoring unit for monitoring the at least one electrical property of the 
organ based on the filtered signal. 

According to still a fiirther aspect of embodiments of the present invention there 
is provided a system for monitoring at least one electrical property of an organ of a 
subject The system comprises a radiofrequency generator for generating an output 
radiofrequency signal and a plurality of electrodes, designed to be connectable to the 
skin of the subject, and configured for transmitting the output radiofrequency signal to 
the organ and sensing an input radiofrequency signal from the organ. The system further 
comprises a monitoring apparatus, eg., the apparatus described herein. 

According to further features in preferred embodiments of the invention 
described below, the dynamically variable filter is adapted in response to a change in a 
physiological condition of the subject The filter is typically a band pass filter 
characterized by a frequency band defined by, e.g., a lower frequency bound and an 
upper frequency bound. 
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According to some embodiments of the present invention the physiological 
condition is a heart rate of die subject 

According to some embodiments of the present invention at least one of a lower 
bound of the frequency band and an upper bound of (he frequency band parameter is a 
linear function of the heart rate. 

According to frirther features in preferred embodiments of the invention 
described below, a lower boimd of the frequency band is about Fl(HR) where 
Fl(HR) = 0,9x(HR/60) Hz and HR is a heart rate of the subject in units of beats per 
minute. 

According to some embodiments of the present invention an upper bound of the 
frequency band is about Fu(HR), where Fu(HR) = 6 + L5x[(HR/60) - 1] Hz. 

According to some embodiments of the present invention the heart rate is 
determined from an ECG signal received from the subject 

According to some embodiments of the present invention the upper frequency 
bound is determined using an iterative process. 

According to some embodiments of the present invention the iterative process is 
based on a comparison between a value of a physiological parameter as extracted from 
the filtered input signal and a value of the physiological parameter as extracted from a 
reference signal. According to some embodiments of the present invention the reference 
signal comprises the ECG signal. According to some embodiments of the present 
invention the physiological parameter is a ventricular ejection time (VET). 

According to some embodiments of the present invention each iteration of the 
iterative process comprises: if the comparison meets a predetermined criterion, then 
updating the upper frequency bound by calculating an average between a low threshold 
for the upper bound and a high threshold for the upper bound. The thresholds can be 
predetermined or they can be set in a previous iteration of the iterative process. 

According to some embodiments of the present invention the physiological 
According to some embodiments of the present invention the iterative process is 
terminated if a value of the VET as extracted from the fdtered input signal is higher than 
a value of the VET as extracted from the reference signal. 

According to some embodiments of the present invention the VET is averaged 
over a plurality of heart beats. 
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According to some embodiments of the present invention the VET is extracted 
&om an average heart beat morphology of the subject. 

Accoidmg to some embodiments of the present mvention an mitial value of the 
upper fiequency bound in the iterative process is a linear function of the heart rate. 
According to some embodiments of the present invention the linear function is FuCHDR). 

According to some embodiments of the present invention the radiofiequency 
signal is filtered using an analog filter. 

According to some embodiments of the present invention at least one quantity is 
calculated using the electrical property. The quantity can be a stroke volume, a cardiac 
output, a brain intra luminal blood volume, a blood flow and the like. 

Accordmg to some embodiments of the present invention the blood flow 
comprises at least one of: an external carotid blood flow rate, an internal carotid blood 
flow rate, an ulnar blood flow rate, a radial blood flow rate, a brachial blood flow rate, a 
common iliac blood flow rate, an external iliac blood flow rate, a posterior tibial blood 
flow rate, an anterior tibial blood flow rate, a peroneal blood flow rate, a lateral plantar 
blood flow rate, a medial plantar blood flow rate and a deep plantar blood flow rate. 

According to some embodiments of the present mvention a phase shift of the 
input radiofiequency signal relative to an output radiofiequency signal transmitted to the 
organ is determined. The phase shift can be used for calculating tiie quantity or 
quantities. 

According to some embodiments of flie present invention the amplitude 
modulation of the input radiofiequency signal is reduced or eUminated so as to provide 
signals of substantially constant envelope. 

According to some embodunents of tiie present invention a phase modulation of 
the input radiofiequency signal is maintained while reducing or eliminatmg the 
amplitude modulation. 

According to some embodiments of tiie present invention tiie input 
radiofiequency signal and fl« output radiofiequency signal are mixed so as to provide a 
mixed radiofirequency signal. According to some embodiments of the present invention 
tiie mixed radiofirequency signal comprises a radiofiequency sum and a radiofiequency 
diffraence. 
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According to some embodiments of the present invention the input 
radio&equency signal is indicative of impedance the organ. According to some 
embodiments of the present invention die input radiofirequency signal is indicative of 
hemodynamic reactance of the organ. 
5 Embodimmts of the present invention successfully address the shortcomings of 

the presently known configurations by providing techniques for processing a signal 
pertaining to one or more electrical property of an organ of a subject 

Unless otherwise defined, all technical and scientific terms used herein have the 
same meaning as commonly understood by one of ordinary skill in the art to which this 

10 invention belongs. Although methods and materials similar or equivalent to those 
described herein can be used in the practice or testing of the present invention, sviitable 
methods and materials are described below. In case of conflict, the patent specification, 
including definitions, will control. In addition, the materials, methods, and examples are 
illustrative only and not intended to be limiting. 

15 Implementation of the method and system of the present invention involves 

performing or completing selected tasks or steps manually, automatically, or a 
combination thereof. Moreover, according to actual instrumentation and equipment of 
preferred embodiments of the method and system of the present invention, several 
selected steps could be implemented by hardware or by software on any operating 

20 system of any firmware or a combination thereof For example, as hardware, selected 
steps of the invention could be implemented as a chip or a circuit. As software, selected 
steps of the invention could be implemented as a plurality of software instructions being 
executed by a computer using any suitable operating system. In any case, selected steps 
of the method and system of the invention could be described as being performed by a 

25 data processor, such as a computing platform for executing a plurality of instructions. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the 
accompanying drawings. With specific reference now to the drawings in detail, it is 
30 stressed that the particulars shown are by way of example and for purposes of illustrative 
discussion of the preferred embodiments of the present invention only, and are presented 
in the cause of providing what is believed to be the most useful and readily understood 
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description of the principles and conceptual aspects of the invention. In this regard, no 
attempt is made to show structural details of the invention in more detail than is 
necessary for a fundamental und^tanding of the invention, the description taken vdth 
the drawings making apparent to those skilled in the art how the several forms of the 
5 invention may be embodied in practice. 
In the drav^gs: 

FIG. 1 is a schematic illustration of a filtering device, according to various 
exemplary embodiments of the present invention; 

FIGs. 2a-b show a representative example of dynamically varying frequency 
10 bounds, employed according to embodiments of tiie present invention; 

FIG. 2c show a representative example of a dynamically varying frequency band, 
employed according to embodiments of the present invention; 

FIG. 3a is a flowchart diagram of a method suitable for processing an input 
signal, according to various exemplary embodiments of the present invention; 
15 FIG. 3b is a flowchart diagram of an iterative process for selecting a frequency 

bound, according to various exemplary embodiments of the present invention; 

FIG. 3c is a schematic illustration of a procedure for extracting ventricular 
ejection time fit>m an EGG signal; 

FIG. 3d is a schematic illustration of a procedure for extracting ventricular 
20 ejection time from a derivative of a filtered signal according to various exemplary 
embodiments of the present invention; 

FIGs. 3e-f are schematic illustrations of beat morphology characterization 
procedures, according to various exemplary embodiments of the present invention; 

FIG. 4 is a schematic illustration of apparatus for monitoring one or more 
25 electriceil properties of an organ of a subject, according to variouis exemplary 
embodiments of the present invention; 

FIG. 5 is a schematic illustration of a signal processing unit, according to various 
exemplary embodiments of the present invention; 

FIG, 6 is a block diagram of electronic circuitry, according to various exemplary 
30 embodiments of the present invention; 
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FIG. 7 is a schematic illustration of a system for monitoring at least one electrical 
property of an organ of a subject, according to a preferred embodiment of the present 
invention; 

FIGs. 8a-e are schematic illustrations showing perspective (Figure Sa), fiont 
5 (Figure 8bX rear (Figure 8c), side (Figure 8d) and top (Figure 8e) views of a sticker 
which can be used for transmitting and sensing the radiofrequency signal, according to 
one embodiment of the present invention; 

FIG. 8f is a schematic illustration of a package of several stickers accor^ng to an 
embodiment of the present invention; 
10 FIGs. 9a-d are schematic illustrations of various electronic circuitries, according 

to exemplary embodiments of the present invention; 

FIGs. lOa-lO-e, lla-lle, 12a-12e, 13a-13e, 14a-14e, I5a-15e, 16a-16gand 17a- 
17g show snapshots of the display of a prototype system, manu&ctured and configured 
according to various exemplary embodiments of the present invention; and 
15 FIGs. 18a-18b, 19a-19b and 20a-20b are plots of cardiac output as calculated 

from a signal filtered according to various exemplary emboctiments of the present 
invention. 

DESCRIPTION OF EXEMPLARY EMBODIMENTS 

20 The present embodiments comprise a method, device apparatus and system 

which can be used for processing a signal. Specifically, but not exclusively, the present 
embodiments can be used for processing a radiofrequency signal sensed from an organ 
of a subject and for monitoring one or more electrical properties of an organ, eg., for the 
purpose of determining one or more quantities which are related to electrical properties. 

25 Thus, for example, exemplary embodiments of the present invention can be used for 
calculating stroke volume, cardiac output, brain intra luminal blood volume and/or blood 
flow. Embodiments of the present invention can also be used for discriminating tissue 
and/or determining at least one of: body cell mass, fat free mass, total body water, 
hematocrit of blood, hydration status and circumference of a body segment. 

30 The principles and operation of a method, device apparatus and system according 

to the present embodiments may be better understood with reference to the drawings and 
accompanying descriptions. 
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Before explaining at least one embodiment of the invention in detail, it is to be 
understood &at the invention is not limited in its application to the details of 
construction and the arrangement of the components set forth in the following 
description or illustrated in the drawings. The invention is capable of other 

5 embodiments or of being practiced or carried out in various ways. Also, it is to be 
understood that the phraseology and terminology employed herein is for the purpose of 
description and should not be regarded as limiting. 

Computer programs implementing the method according to embodiments of the 
present invention can commonly be distributed to users on a distribution medium such 

10 as, but not limited to, a floppy disk, CD-ROM and flash memory cards. From the 
distribution medium, the computer programs can be copied to a hard disk or a similar 
intermediate storage medium. The computer programs can be run by loading the 
computer instructions either from their distribution medium or their intermediate storage 
medium into the execution memory of the computer, configuring the computer to act in 

IS accordance with the method of this invention. All these operations are well-known to 
those skilled in the art of computer systems. 

A typical system for monitoring electrical properties of a body section, such as a 
bioimpedance system, includes a tetrapolar array of circumferential band electrodes 
cormected to the subject at the base of the neck and surrounding the circumference of the 

20 lower chest, at the level of die xiphoid process. When a constant magnitude alternating 
current flows through the upper cervical and lower thoracic band electrodes, a voltage, 
proportional to the thoracic electrical impedance (or reciprocally proportional to the 
admittance), is measured between the irmer cervical and thoracic band electrodes. The 
portion of the cardiac synchronous impedance change, temporally concordant with the 

25 stroke volume, is ascribed solely and uniquely to volume changes of the aorta during 
expansion and contraction over the heart cycle. A typical printed circuit board of such 
system comprises one or more band pass filters, a half-wave rectification circuit and one 
or more low pass filters. 

The present Inventor discovered techniques for reducing the noise associated 

30 with conventional systems. As demonstrated in the Examples section that follows, the 
present Inventor succeeded in reducing noise introduced due to patient agitation or other 
physiological phenomena like breathing. The present Inventor discovered techniques for 



wo 2009/022330 



PCT/IL2008/001105 



11 

separating and differentiating between a cardiovascular bioreactance signal and a 
respiratory bioreactance signal, where the latter is typically much larger than the former. 

The present Inventor has realized that tiie noise level is proportional to the 
bandwidth of the band pass filter and that a considerable portion of the noise passes the 

S band pass filter hence being folded into the half-wave rectification circuit 

The present Inventor also discovered techniques for reducing or eliminating AM 
noise hence significantly improving the ability to provide accurate measurement 

Referring now to the drawings. Figure 1 illustrates a filtering device 10, 
according to various exemplary embodiments of the present invention. Device 10 

10 comprises a first input unit 12 which receives an input signal 16 pertaining to one or 
more electrical properties of an organ of a subject. For example, signal 16 can relate to 
the hemodynamic reactance of the organ. 

As used herein, "hemodynamic reactance" refers to the ims^inary part of the 
impedance. Techniques for extra^cting the imaginary part fix>m the total impedance are 

15 known in the art Typically, such extraction is perforaied at hardware level but the use 
of algorithm at a software level is not excluded firom the scope of the present invention. 
Signal 16 can be provided, for example, by processing a radiofi-equency signal sensed 
firom the organ, as fiirtiier detailed heieinunder. 

In various exemplary embodiments of the invention device 10 fijrther comprises 

20 a second input unit 14 which receives data 18 pertaining to a physiological condition of 
the subject The physiological condition is preferably, but not obligatorily, the heart rate 
of the subject, and the data pertaining to the physiological condition can be analog data 
or digital data, as desired. While the embodiments below are described with a particular 
emphasis to physiological condition which is a heart mte, it is to be understood tliat 

25 more detailed reference to the heart rate is not to be interpreted as limiting the scope of 
the invention in any way. For example, in exemplary embodiments of the present 
invention the physiological condition is a ventilation rate of the subject, a repetition rate 
of a peirticular muscle unit and/or one or more characteristics of an action potential 
sensed electromyography. 

30 Device 10 further comprises a filtering unit 20 which filters the input signal 16 to 

provide a filtered signal 22. In various exemplary embodiments of the invention the 
filtering is according to a firequency band which is dynamically adapted in response to a 
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change in fhe physiological condition of the subject. It was found by the Inventor of the 
present invention that the dynamical adaptation of the frequency band to the 
physiological condition of the subject can significantly reduce the influence of unrelated 
signals on the measured or monitoring of electrical properties of the body section. 
5 The adaptation of the frequency band to the physiological condition can be 

according to any adaptation scheme known in the art. For example, one or more 
parameters of the fiiequency band (e.g., lower bound, upper bound, bandwidth, central 
frequency) can be a linear fimction of a parameter charact^zing the physiological 
condition. Such parameter can be, for example, the number of heart beats per minute. 

10 A representative example of a dynamically varying frequency bounds, employed 

according to some embodiments of the present invention by unit 20, is illustrated in 
Figures 2a-b. Shown in Figures 2a-b is the frinctional dependence of the frequency 
bounds (upper bound in Figure 2a and lower bound in Figure 2b) on the heart rate of the 
subject. As shown in Figure 2a, the upper bound of the frequency band varies linearly 

15 such that at a heart rate of about 60 beats per minute (bpm) the upper bound is about 6 
Hz, and at a heart rate of about 180 bpm the upper bound is about 9 Hz. As shown in 
Figure 2b, the lower bound of the frequency band varies linearly such that at a heart rate 
of about 60 the lower bound is about 0.9 Hz bpm and at a heart rate of about 180 bpm 
the lower bound is about 2.7 Hz. 

20 As used herein fhe terai "about" or "approximately" refers to ± 10 %. 

In some embodiments of the present invention the upper bound approximately 
equals the function Fu(HR) defined as Fu(HR) = 6 + 1.5 x [(HR/60) - l]Hz, where HR is 
the heart rate of the subject in units of bpm. In some embodiments, the upper bound 
equals Fu(HR) at all times, while in other embodiments, the upper bound is set using an 

25 iterative process. 

In some embodiments of the present invention the lower bound approximately 
equals the fimction Fl(HR) defmed as Fl(HR) == 0.9 x (HR/60) Hz. In some 
embodiments, the lower bound equals Fl(HR) at all times while in other embodiments 
the lower bound is set by an iterative process. 

30 Representative examples of iterative process suitable for some embodiments of 

the present invention are provided hereinunder. 
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A dynamically varying band pass filter (BPF) characterized by a dynamically 
varying upper frequency bound and a dynamically varying lower frequency boimd, 
according to some embodiments of the present invention is illustrated in Hgure 2c. As 
shown, each heart rate is associated with a frequency band defined by a lower bound and 

5 an upper bound. For example, for a heart rate of 60 bpm. Figure 2c depicts a BPF in 
which the lower bound is about 0.9 Hz and the upper bound is about 6 Hz. 

It is to be understood that the values presented above and the functional relations 
illustrated in Figures 2a-b are exemplary embodiments and should not be considered as 
limiting the scope of the present invention in any way« In other exemplary 

10 embodiments, the fimctional relations between the frequency band and the physiological 
condition can have different slopes and/or offsets, or they can be non-linear. 

Figiire 3a is a flowchart diagram of a method suitable for processing an input 
signal, according to various exemplary embodiments of the present invention. The 
method can be executed by activating filtering device^ 10 or using any other filt^ing 

15 device supplemented by appropriate circuitry. Selected steps of the method can be 
embodied in many forms. For example, the selected steps can be embodied in on a 
tangible medium such as a computer for performing the selected steps. The selected 
steps can be embodied on a computer readable medium, comprising computer readable 
instructions for carrying out the selected steps. The selected steps can also be embodied 

20 in electronic device having digital computer capabilities arranged to run the computer 
program on the tangible medium or execute the instruction on a computer readable 
medium. 

The method begins at step 30 and continues to step 31 in which the physiological 
condition of the subject is determined. The physiological condition can be, as stated, a 

25 heart rate and it can be determined using any procedure known in the art, such as, but 
not limited to, analysis of an electrocardiogram (ECG) signal or the like. The method 
continues to step 32 in which a frequency band is selected based on the physiological 
condition of the subject, and proceeds to step 33 in v^iiich the input signal is filtered 
according to frequency band. In various exemplary embodiments of the invention the 

30 method loops back to step 31 so as to dynamically adapt the frequency band in response 
to a change in the physiolo^cal condition. 
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Hie selection of frequency band can be according to any adaptation sch^e, 
including, without limitation, the use of one or more linear functions (e.g.^ the functions 
Fu and Fl) as further detailed hereinabove, and/or an it^tive process as furdier detailed 
hminbelow. 
5 The method ends at step 34. 

Following is a description of an iterative process for determining the fiequency 
band of tiie band pass filter which filters the input signal according to some 
embodiments of the present inventioxu The iterative process can, in some embodiments, 
based a comparison between a value of a physiological parameter as extracted or 
10 calculated fix)m the filtered input signal and a value of the same physiological parameter 
as extracted or calculated firom a reference signal, for example, an BCG signal. 

Hie term physiologjlcal parameter refers to any variable paramet^ which is 
measurable or calculable and is representative of a physiological activity, particularly, 
but not necessarily, activity of the heart In various exemplary embodiments of the 
15 iiivention the physiological parameter is other tiian the heart rate per se. The 
physiolo^cal parameter can be a time-related parameter, amplitude-related parameters 
or combination thereof. 

Typically, the filter signal and the reference signal are expressed in terms of 
amplitude as a function of the time. Thus, time-related parameters are typically 
20 calculated using abscissa values of the signals and amplitude-related parameters are is 
typically calculated using ordinate values of the signals. 

Representative of time-related physiological parameters suitable for the present 
embodiments include, without limitation, systolic time, diastolic time, pre-ejection 
period and ejection time. A representative example of amplitude-related physiological 
25 parameter suitable for the present embodiments includes, without limitation, cardiac 
contractility, maximal amplitude above zero during a single beat, maximal peak-to-peak 
amplitude during a single beat, and RMS level during a single beat Also contemplated 
are various slopes parameters, such as, but not limited to, the average slope between two 
points over the signal. 

30 In various exemplary embodiments of the invention the physiological parameter 

is a ventricular ejection time (VET). 
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While the embodiments below are described with a particular emphasis to VET 
as the physiological parameter, it is to be understood that more detailed refer^ce to 
VET is not to be inteipreted as limiting the scope of the invention in any way. 

The present inventors discovered that a significant amount of the biological 
S information for a particular subject can be obtained from a frequency range between 
Fl(HR) and 5.5 Hz, where HR is the heart rate of the subject It was further discovered 
by the present inventors that for some medical conditions some of the information can 
reside between 5.5 Hz and Fu(HR). 

The advantage of the comparison between two different techniques for extracting 
10 or calculating the same physiological parameter, is that it allows to substantially 
optimize the upper frequency bound of the of the band pass filter. In various exemplary 
embodiments of the invention in each iteration of the itemtive process, the comparison is 
repeated. If the comparison meets a predetermined criterion, the upper frequency bound 
is updated by calculating an avenge between a low threshold for the upper bound and a 
15 high threshold for the upper bound. The lower frequency bound can be a constant 
bound, e.g.^ a constant frequency \^ich is from about 0.9 Hz to about 2.7 Hz), or it can 
be dynamic, e.g., Fl(HR), HR being the heart rate of the subject before or during the 
respective iteration. 

Ihe low and high thresholds for the upper bound can be set in more than one 
20 way. In some embodiments, the low and high thresholds are predetermined (namely 
they determined a priori before the iterative process), in some embodiments, the 
thresholds are set in a previous iteration of iterative process, in some embodiments one 
of the thresholds is predetermined and the other threshold is set in a previous iteration of 
iterative process. In any event, the first iteration is based on two thresholds which are 
25 determined a priori before the iterative process. It was found by the inventors of the 
present invention that, at least initially (/.e., at the first iteration), the first threshold can 
be about Fu(40), which in various exemplary embodiments of the invention is about 
5.5 Hz, and the second threshold can be the calculated value of Fu(HR), HR being the 
heart rate of the subject before or during the respective iteration. 
30 The predetermined criterion used during the iterations can be, for example, that 

the results of the two calculations are similar (e.g.^ within about 40 % or 30 % or 25 % 
of each other). Ihe predetermined criterion can also relate to the direction of difference 
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between fhe two calculations. Broadly, for time-Telated parameters, the upper bound is 
updated if the value of the parameter as calculated based on the reference signal is 
higher than value of the parameter as calculated based on the filtered signal, and for 
amplitude-related parameters the upper bound is updated if tilie value of parameter as 
5 calculated based on the reference signal is lower than the value of the paramet^ as 
calculated based on the filtered signal. For slope-related parameters, the upper bound is 
typically updated if the value of the parameter as calculated based on tfie reference 
signal is higher than the value of the parameter as calculated based on the filtered signal. 
A Boolean combination between the above criteria can also be used as a 

10 criterion. For example, an AND Boolean combination can be employed in which case 
the upper fi:equency bound can be updated if the results of the two calculations are 
similar and the calculation according to the filtered signal indicates an abnormal 
physiological condition while the calculation according to the reference signal indicates 
a normal physiological condition. 

15 Figure 3b is a flowchart diagram of an itemtive process for selecting the upper 

firequency bound, according to various exemplary embodiments of the present invention. 
The description is for a physiological parameter which is VET, but, as stated, it is not 
intended to limit the scope of the present invention to this ^e of physiological 
parameter. 

20 The iterative process begins at 60 and continues to 61 in which the upper 

frequency bound is set to the value of Fu(HR), HR bemg the heart rate of the subject 
before the initiation of the iterative process. The heart rate can be inputted or it can be 
determined by the process, eg., firom the ECG signal. The process continues to 62 in 
which initial values are assigned to two frequency thresholds. The frequency thresholds 

25 are denoted in Figure 3b by Tl and T2. In various exemplary embodiments of the 
invention the initial value of Tl is Fu(40) and the initial value of T2 is the initial upper 
frequency bound. 

The process continues to 63 and 64 in which VET is calculated separately from 
the filtered input signal (63) and from a reference signal (64), e.g,, ECG. The VET as 
30 calculated from the input signal is denoted in Figure 3b VETl and the VET as calculated 
from the ECG signal is denoted in Figure 3b VET2. In the first iteration, the filtered 
signal from which VETl is calculated is preferably obtained by filtering the input signal 
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using a band pass filter defined between the upper and lower fiiequency bounds as 
initially set at 61. 

The process continues to decision 65 in which VETl and VET2 are compared. If 
VETl is higher than VET2, the iterative process continues to 73 where it is terminated^ 
5 If VETl is not higher than VET2, the process continues to 66 in wiiich the upper bound 
is updated to the average AVE of Tl and T2. In various exemplary embodiments of the 
invention .AVE(T1, T2) is the arithmetic mean of Tl and T2 (ie., AVE(T1, 
T2) = (Tl+T2)/2), tWs need not necessarily be the case, since, in some embodiments, a 
different averaging scheme (e.g:, a weighted average, a geometric mean, harmonic 
10 mean, RMS, etc.) can be employed. 

From 66 the process continues to 67 in which VETl is recalculated but from a 
signal which is filtered using an updated band pass filter. The lower bound of the 
updated band pass filter can be the initial lower bound or it can be updated based on the 
heart rate of the subject immediately before the filtration of the input signal (e.g., 
15 according to Fl described above). The upper bound of the updated band pass filter is 
preferably the updated upper bound. Optionally, the process continues to 68 in which 
the VET2 is also recalculated from the reference signal. Alternatively, the value of 
VET2 from 64 can be used. 

The process then continues to decision 69 in wiiich the process, determines 
20 whether or not a predetermined termination criterion is met. If the piedetermined 
termination oriterion is met, the iterative process continues to 73 where it is terminated, 
otherwise the process continues to decision 70 in which VETl and VBT2 are compared. 
If the deviation between VETl and VET2 is lower than a predetermined threshold A and 
VETl is lower than VET2, the process continues to 71 at vsiiich the value of the upper 
25 bound is assigned to T2, otherwise the process continues to 72 at which the value of the 
upper bound is assigned to TL 

From 71 and 72 the process loops back to 66 at which an additional iteration 

begins. 

The threshold A employed at decision 70 is typically expressed as a fraction of 
30 VET2 or a fraction of (VET2-VET1)A/^ET2. In various exemplary embodiments of flie 
invention A = p* VET2, where p < 0.5, for example, p = 0.4 or p = 0,3 or p = 0.25. 



wo 2009/022330 



PCT/IL2008/001105 



18 

The termination criterion employed at decision 69 can be for example, a maximal 
number of iterations. In this embodiment, the process counts the number of iterations 
and compares them to a predetemiined iteration number threshold. If the number of 
iterations exceeds iteration number Areshold the process determines that the termination 
5 condition is met and continues to end 73. Typically, but not necessarily, the value of the 
predetermined iteration nwnber threshold is from about 3 iterations to about 10 
iterations, e.g,, 5 or 6 iterations. 

The termination criterion can also include a comparison of the value of VETl to 
a predetermined VET threshold Tvet, which may be absolute, subject-specific or relative 

ID to VET2. It this embodiment, the process compares the value of VETl to Tvet- If 
VETl is higher than Tvet the process determines that the termination condition is met 
and continues to end 73. Typically, but not necessarily, a relative value of Tvct is 
employed. For example, in various exemplary embodiments of the invention 
TvET = p* VET2 where p < 0.5, for example, p = 0.4 or p = 0.3 or p = 0.25. 

15 In any of the above embodiment, the VET (either VETl or VET2) can be 

calculated over a single heart beat, or, more preferably, it can be averaged over two or 
more heart beats. The advantage of using an avera^g procedure rather than using a 
single beat is that it attenuates random disturbances v^iich may be present in the signal 
during the iterative process. Nevertheless, extmction of the VET from a single beat is 

20 not excluded firom the scope of the present invention. 

For a single heart beat, the calcidation of VET can be performed by 
characterizing the morphology of the beat, identifying two or more identifiable points on 
the beat and measuring the time between the identified points. 

Figure 3c illustmtes a procedure for extracting VET2 when the reference signal 

25 is an ECG signal. Shown in Figure 3c is a typical morphology of a single beat of an 
ECG signal as a function of the time. VET2 can be defined as the time period 
(difference between the abscissa values) between the R peak and the T peak of the ECG 
signal. 

When the filtered signal is hemodynamic reactance, the value of VETl is 
30 preferably extracted fi-om the first derivative of the filtered signal. The procedure is 
illustrated in Figure 3d, which illustrates a typical morphology of a single beat of the 
hemodynamic reactance N and its first derivative dW, as a function of the time. As 
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shown, <W has two zeroes Oi and O2 over the beat, with a pomt of local mayimum Mi 
between tiie zeroes and a point of local minimum M2 after the second zero. In some 
embodiments of the present invention VETl is defined as the time period (difference 
between the abscissa values) between the first zero Oi and the first minlmiiTu M2 after 
5 the second zim> O2. 

The beat morphology of the filtered signal can be characterized by identifying 
identifiable points on the reference signal and using the time associated with these points 
(abscissa values) for defining anchor points on the filtered signal. 

Two exemplary beat morphology chamcterization procedure of the filtered signal 
10 for embodiments in which the ir^>ut signal is the hemod3mamic reactance and the 
reference signal is the ECG signal are illustrated in Figures 3e-f. 

In Figure 3e, a single beat of d^ (first derivative of the hemodynamic reactance 
N) is defined between two anchor endpoints: a first (left) endpoint has the abscissa value 
of the Q peak of the ECG signal, and a second (right) endpoint has the abscissa value of 
15 the R peak of the next ECG signal. In other words, a single beat of diV^ has a width 
which equals the following sum of five successive intervals over the ECG: 
QR + RS + ST + TQ + QR. 

In Figure 3£ a single beat of d^ is defined using three anchor points: the two 
anchor endpoints as described in Figure Be and an intermediate anchor point which has 
20 the abscissa value of the global maximum A of the hemodynamic reactance N between 
the two endpoints. 

When the VET is averaged, the calculation can be done in more than one way. 
In some embodiments, the same morphology characterization is employed for a 
plurality of beats over a predetermined interval of the respective signal, so as to provide 
25 one local VET for each beat. The VET can be defined as the average of all local VETs. 
Any averaging procedure can be employed, include, without limitation, arithmetic mean, 
weighted average, geometric mean, harmonic mean, RMS and the like. 

In some embodiments, the morphologies of the beats are averaged over an 
ensemble of beats in the filtered signal, to provide an average beat morphology, and the 
30 VET is be determined by identifying two or more identifiable points on the average beat 
morphology and measuring the time between the identified points. For example, for 
each beat in the ensemble, the morphology can be characterized as described above and 
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all the morphologies can be averaged, e.g:, pomt by point. The morpholo^es can also 
be averaged in segments. This embodiments is particularly useful when a single beat is 
defined using more than two anchor points in which case each segment over the beats 
(between two successive anchor points) can be averaged, e.g.^ point-by-point The 
average beat morphology can then be obtained by stitching the averaged segments. For 
example, in the embodiment illustrated in Figure 3f the beat is defined using two 
intermediate anchor points and an intermediate point In this embodiment, each beat has 
a left segment (firom the abscissa value of Q to the abscissa value of A) and a right 
segment (fix)m the abscissa value of A to the abscissa value of R). Hie left segments of 
all beats in the ensemble can be averaged to provide a left segment average, and the right 
segments of all beats in the ensemble can be averaged to provide a right segment 
average. The average beat morphology can be obtained by stitching the left segment 
average to the right segment average. 

In varioiis exemplary embodiments of the invention the time scale of the beats in 
the ensemble is adjusted so as to fit all the beats in the ensemble to a single time scale. 

The result of this average is a single beat morphology from which VETl can be 
extracted as described above. For example, when the first derivative of the 
hemodynamic reactance is used, the average morphology typically has the shape 
illustrated in Figure 3d. VETl can then be extracted as the time period between 0| and 
02- The first derivative of the signal can be calculated before or after averaging. When 
the first derivative is calculated before the averaging, the averaging procedure described 
above is performed with respect to the derivative of the signal. When the first derivative 
is calculated after the averaging, the averaging procedure described above is performed 
with respect to the signal and the obtained average is then differentiated. Calculation of 
the first derivative after averaging is preferred fi*om the standpoint of noise reduction. 

The predetermined time period over which an averaging procedure is performed 
to extract the VET typically extends over about 10 heart beats. 

A particular advantage of the device and method of the present embodiments is 
that they can be implemented in many systems designed for measuring or monitoring 
electrical properties of body sections, thereby improving their performance, e.g., by 
increasing their signal to noise ratio at least for situations in which the amount of noise 
is high. Representative examples of such systems include, without limitation, a system 
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for monitoring blood flow, cardiac output and/or stroke volume, which can be similar to 
or based on the systems disclosed in U.S. Published Application No. 2006020033 and 
International Patent Publication No. WO2006/087696, the contents of which are hereby 
incorporated by reference; a system for predicting body cell mass, fat free mass and/or 

5 total body water of a subject, which can be similar to or based on the system disclosed in 
U.S. Patent No. 5,615,689, the contents of which are hereby incorporated by reference; a 
system for determining hematocrit of blood in a body part of a subject, which can be 
similar to or based on the system disclosed in U.S. Patent No. 5,642,734, the contents of 
vMch are hereby incorporated by reference; a system for monitoring hydration status of 

10 a subject, which can be similar to or based on the system disclosed in U.S. Published 
Application No, 20030120170, the contents of which are hereby incorporated by 
reference; a system for discriminating tissue, which can be similar to or based on the 
system disclosed in U.S. Published Application No. 2006008S048» the contents of \^ch 
are hereby incorporated by reference; and a system for calculating the circumference of 

15 a body segment which can be similar to or based on the system disclosed in U.S. 
Published Application No. 20060122540, tiie contents of which are hereby incorporated 
by reference. 

Reference is now made to Figure 4 which is a schematic illustration of apparatus 
40 for monitoring one or more electriced properties of an organ of a subject 121, 

20 according to various exemplary embodiments of the present invention. 

Apparatus 40 comprises an input unit 42 for receiving an input radiofrequency 
signal sensed from the organ. The input radiofrequency signal typically comprises a 
radiofrequency signal related to the electrical properties of the organ {e.g., bioimpedance 
which may generally relate to the impedance and/or hemodynamic reactance of the 

25 organ). The signal is sensed from one or more sensing locations 48 on the organ of 
subject 121 and is originated from an output radiofrequency signal 124 generated by a 
radiofrequency generator 122. The input radiofrequency signal, however, can include 
one or more noise components, which may be introduced into the signal due to various 
reasons, e.g,^ subject agitation or breathing. In various exemplary embodiments of the 

30 invention apparatus 40 is capable of reducing or eliminating these noise components. 

Apparatus 40 fiuther comprises a signal processing unit 44 which processes the 
input radiofrequency signal. The processing may include, for example, mixing. 
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demodulation, determination of phase shift, analog filtering, sampling and any 
combination thereof. Signal processing unit 44 may or may not be in communication 
with radiofrequency generator 122, as desired. A representative example of signal 
processing unit 44 is provided hereinunder with reference to Figure S. 
5 Apparatus 40 further comprises a filtering unit 46 vAnch filters the processed 

input signal. Unit 46 preferably perfomis the filtration operation in the frequency 
domain. Tlius, in various exemplary embodiments of the invention, a series of samples 
of the processed radiofrequency signal is transformed, e.g.^ by a Fast Fourier Transform 
(FFT), to provide a spectral decomposition of the signal in tiie firequency domain. The 
10 transformation to the frequency domain can be done by a data processor. Algorithms for 
performing such transformations are known to those skilled in the art of signal 
processing. 

The obtained spectral decomposition of the signal is filtered by unit 46 which 
typically eliminates one or more of the frequencies in the spectrum, depending on the 

15 upper and lower frequency bounds of the filter employed by unit 46. Unit 46 preferably 
employs a dynamically variable filter. For example, unit 46 can comprise filtering 
device 10 as described above. 

Once filtered, the signal is transmitted to a monitoring unit 52 which monitors 
the electrical property or properties of the organ based on filtered signal. Unit 52 can 

20 monitor the electrical property by recording it and/or transmitting it to an external 
device, such as a display device and/or a computer. The dynamically variable filter can 
be adapted in response to a change in the physiological condition of the subject, as 
fiuther detailed hereinabove. 

Apparatus 40 is optionally and preferably designed for determining a phase shift 

25 of signal 126 relative to signal 124. This can be done using a phase shift 

determinator 50 (not shown, see Figure 5) which can operate according to any known 
technique for determining a phase shift. The phase shift can be determined for any 
frequency component of the spectrum of radiofrequency signals received from the 
organ. For example, in one embodiment, the phase shift is determined from the base 

30 frequency component, in another embodiment the phase shift is determined from the 
second frequency component, and so on. Alternatively the phase shift can be 
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detennined using several frequency components, e.^., using an appropriate averaging 
algorithnt 

It was discovered by the Inventor of the present invention that the phase shift of 
the ii4>ut signal, as received from the organ, relative to the output signal as generated by 
5 generator 122, is indicative of tiie blood flow in the organ. Thus, according to some 
embodiments of the invention the blood flow is detennined using the phase shift (A^). 

The advantage of using for detennining the blood flow is that the relation 
between the blood flow and Ag> depends on fewer measurement-dependent quantities as 
compared to prior art determination techniques in which the impedance is used. 

10 Specifically, it was found by the Inventor of the preset invention that there is a linear 
relationship between Aip and the blood flow, with a proportion coefficient comprising 
the systolic ejection time, T. For example, the stroke volume SV can be calculated using 
the relation SV = const xT^A^?, and the cardiac output CO can be calculated using the 
relation CO = const. x^xA^xHR, where HR is the heart rate of the subject (e,g., in units 

15 of beats per minute), and "const." a constant which can be found, for example, using a 
calibration curve. As will be appreciated by one ordinarily skilled in the art, the absence 
of i and Zo from the formulae for SV and CO significantly reduces the uncertainty in the 
obtained values because there is no entanglement between the obtained values and ertors 
associated with the measurement of L and Zo. 

20 In various exemplary embodiments of the invention apparatus 40 comprises a 

data processor 142, configured for calculating at least one quantity using the filtered 
signal. Data processor 142 can also be employed by unit 46 for performing the 
transformation to the frequency domain and/or eliminating the firequency components 
according to the dynamically variable fi^uency bounds. 

25 Many quantities may be calculated by data processor 142. For example, in 

various Gxemplary embodiments of the invention processor 142 calculates blood-volume 
related quantities, such as, but not limited to, a stroke volume, a cardiac output and a 
brain intra luminal blood volume. In the embodiments in which data processor 142 is 
employed, monitoring unit 46 can monitor the quantity calculated by processor 142. 

30 When apparatus 40 codetermines the phase shift, processor 142 can calculate the 
quantity based on the phase shift. 
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Reference is now made to Figure S which schematically illustrates signal 
processing unit 44, according to various exemplary embodiments of tfie present 
inventioiL Unit 44 preferably comprises a mixer 128, electrically conununicating with 
generator 122, for mixing signal 124 and signal 126, so as to provide a mixed 
5 radiofirequency signal. Signals 124 and 126 may be inputted into mixer 128 through 
more than one channel, depending on optional analog processing procedures 
amplification) which may be performed prior to the mixing. 

Mixer 128 may be any known radio&equency mixer, such as, but not limited to, 
double-balanced radiofirequency mixer and unbalanced radiofrequency mixer. 
10 According to a preferred embodiment of the present invention, the mixed radiofirequency 
signal is composed of a plurality of radiofrequency signals, which may be, in one 
embodiment, a radiofrequency sum and a radiofrequency difference. A sum and a 
difference may be achieved, e.g., by selecting mixer 128 such that signal 124 and signal 
126 are multiplied thereby. Since a multiplication between two frequencies is equivalent 
15 to a frequency sum and a frequency difference, mixer 128 outputs a signal which is 
composed of the desired radiofrequency sum and radiofrequency difference. 

The advantage in the production of a radiofrequency sum and a radiofrequency 
difference is that wdiereas the radiofrequency sum includes both the signal, v\^ch is 
indicative of the electrical property, and a considerable amount of electrical noise, the 
20 radiofrequency difference is approximately noise-free. 

It was found by the present Inventor that this technique is suitable for minimizing 
the electrical noise even when the effect of interest is smaller than the measured quantity 
by 2-4 orders of magnitude. 

According to various exemplary embodiments of the present invention unit 44 
25 further comprises a phase shift determinator 50 for determining the phase shift of the 
input signal relative to the output signal. Phase shift determinator 50 can determine the 
phase shift according to any technique known in the art. For example, the phase shift 
can be determined from the radiofrequency difference outputted from mixer 128. 

According to a preferred embodiment of the present invention processing unit 44 
30 further comprises electronic circuitry 132, which filters out a portion of the signal such 
that a remaining portion of the signal is characterized by a substantially increased signal- 
to-noise ratio. 
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Circuitry 132 is better illustrated in Figure 6. According to an embodiment of 
the present invention circuitry 132 comprises a low pass filter 134 to filter out the high 
frequency content of the signal. Low pass filter 134 is particularly usefiil in the 
embodiment in which mixer 128 outputs a sum and a difiEerence, in which case low pass 
5 filter 134 filters out the ladiofiiequency sum and leaves the approximately noise-fi:ee 
radiofi:equency difference. Low pass filter 134 may be designed and constructed in 
accordance with the radiofirequency difference of a particular system which employs 
apparatus 40. A judicious design of filter 134 substantially reduces the noise content of 
the remaining portion. In a conventional bioimpedance system, for example, a 
10 substantial amount of the noise of the received signal is folded into the remaining signal, 
which is thus characterized by a bandwidth of about 2 kilohertz. It has been found by 
the inventor of the present invention that by including ou^ut radiofi:equency signal 124 
and by mixing it with input radiofirequency signal 126, the noise in the resulting signal is 
characterized by a bandwidth that is at least one order of magnitude below the noise 
15 bandwidth of conventional systems. 

In various exemplary embodiments of the invention mixer 128 and circuitry 132 
are designed and constructed for increasing the signal-to-noise ratio by at least 20 dB, 
more preferably by 25 dB, most preferably by 30 dB. 

Circuitry 132 prefembly comprises an analog amplification circuit 136 for 
20 amplifying the remaining portion of the signal. The construction and design of analog 
amplification circuit 136 is not limited, provided curcuit 136 is capable of amplifying the 
signal. A non limiting example of amplification circuit 136 is further detailed herein 
below in the Examples section that follows. 

According to a preferred embodiment of the present invention circuitry 132 
25 fiirther comprises a digitizer 138 for digitizing the signal. The digitization of the signal 
is useful for furdier digital processing of the digitized signal, e.g., by a microprocessor. 

Optionally, circuitry comprises a differentiator 140 (either a digital differentiator 
or an analog differentiator) for performing at least one time-differentiation of the 
measured impedance to obtain a respective derivative (e,g;, a first derivative, a second 
30 derivative, etc.) of the electrical property. Differentiator 140 may comprise any known 
electronic functionality (e.g., a chip) that is capable of performing analog or digital 
differentiation. Time-derivatives are useful, for example, when the electrical property is 



wo 2009/022330 PCT/IL2008/001105 

26 

bioimpedance and the apparatus is employed in a system for measuring stroke volume or 
cardiac output, as further detailed hereinafter. 

According to a preferred embodiment of fhe present invention signal processing 
unit 44 comprises an envelope elimination unit 135 which reduces or, more preferably, 
5 eliminates amplitude modulation of signal 126. Optionally and preferably, unit 135 
maintains the phase modulation of signal 126. The input to envelope elimination unit 
135 typically carries a substantial amount of AM noise, which can be described, wittiout 
limitation as a signal vze = v(/)cos(^+^/)), which contains both phase and amplitude 
modulation. According to a preferred embodiment of the present invention unit 135^ 

10 generates a signal having a substantial constant envelope, e.g.^ V26' = voCOs(a7/4-^r)), 
where vo is substantially a constant. The output of unit 135 thus represents the phase (or 
frequency) modulation of signal 126. Unit 135 can employ, for example, a limiter 
amplifier which amplifies signal 126 and limits its amplitude such that the amplitude 
modulation is removed. The advantage of the removal of the amplitude modulation is 

15 that it allows jsl better determination of the phase shift Atp between the input and output 
signals, as further detailed hereinabove. 

Reference is now made to Figure 7, which is a schematic illustmtion of system 
120 for monitoring at least one electrical property of an organ of a subject, according to 
a preferred embodunent of the present invention. System 120 preferably comprises a 

20 radiofrequency genemtor 122, for generating an output radiofrequency signal. 
Generator 122 may be embodied as any radiofrequency generator, such as, but not 
limited to, radiofrequency generator 112 of system 110. System 120 further comprises a 
plurality of electrodes 125, which are connected to the skin of subject 121. Electrodes 
125 transmit output radiofrequency signal 124, genemted by generator 122 and sense 

25 input radiofrequency signal 126 originated from the organ of subject 121. 

System 120 preferably comprises any of the components of apparatus 40 
described above. According to a preferred embodiment of the present invention system 
120 further comprises a detector 129 for detecting a voltage drop on a portion of the 
body of subject 121 defined by fhe positions of electrodes 125. In response to the 

30 detected voltage, detector 129 preferably generates a signal which is indicative of 
impedance of the respective portion of the body. In this embodiment, the stroke volume 
can be calculated using (cU7^0niax» as further detailed hereinabove. Knowing the stroke 
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volume, the cardiac ou^ut is calculated by multiplying the stroke volume by the heart 
rate of the subject More preferably, detector 129 generates a signal which is indicative 
of a hemodynamic reactance, X. 

The blood flow determination provided by system 120 may be used both for 

5 diagnostic and for treatment. Hence, according to a preferred embodiment of the present 
invention, system 120 may further comprise a pacemaker 144. In this embodiment, the 
data processor (not shown, see Figure 4) is preferably programmed to electronically 
control pacemaker 144 in accordance with the calculated quantity. For example, in one 
embodiment, the data processor calculates the cardiac output and sends signals to 

10 pacemaker 144 which controls, substantially in real-time, the heart rate of subject 121, 
so as to improve the cardiac output. 

Additionally or alternatively, system 120 may also comprise a cardiac assist 
device 148, preferably constmcted and design for increasing the cardiac output. Cardiac 
assist devices are known in the art and typically comprise a reinforcing member which 

15 restricts an expansion of a portion of the heart tissue, so that the cardiac output is 
increased. In this embodiment, the data processor is preferably programmed to 
electronically control device 148 in accordance with the calculated cardiac oulput, so 
that both the determination and the improvement of the cardiac output are automatically 
performed by system 120. 

20 According to a preferred embodiment of the present invention system 120 

comprises a drug administrating device 146. Device 146 serves for administrating drugs 
to subject 121. In this embodiment, the data processor is preferably programmed to 
electronically control device 146, in accordance with the value of the calculated 
quantity. For example, if the calculated quantity is the brain intra luminal blood volume, 

25 then, depending on the value of the blood volume, the data processor sends signal to 
device 146 and thereby controls the amount and/or type of medications administered to 
subject 121, 

It is to be understood that any number of electrodes of system 125 or connection 
configurations of electrodes 125 to subject 121 are not excluded from the present 
30 invention. Any type of electrode, in any combination, may be used, for measuring blood 
flow in any artery of the body, such as, but not limited to, the external carotid artery, the 
internal carotid artery, the ulnar artery, the radial artery, the brachial artery, the common 
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iliac artery, the external iliac artery, the posterior tibial artery, the anterior tibial artery, 
the peroneal artery, the lateral plantar artery, the medial plantar artery and the deep 
plantar artery. 

When system 120 is used together with other systems it is desired to minimize 
5 the area occupied by electrodes 125 so as not to interfere the operation of the other 
systems. For example, in intensive care units, the subjects are oftentimes connected to 
ECG leads, arterial line, central venous line, brain stem evoked response equipment, 
chest tubes, GI tube, intravenous and the like. 

Figures 8a-e are schematic illustrations showing perspective (Figure 8a), fix>nt 

10 (Figure 8b), rear (Figure 8c), side (Figure 8d) and top (Figure 8e) views of a sticker 141 
which can be used for transmitting and sensing the radiofiequency signal, according to 
one embodiment of the present invention. The sticker comprises electrical contacts 145 
being as fixed and predetermined distance therebetween, tiius reducing any the efiTect of 
variable inter-electrode distance on the measurement. Two such contacts 145a and 145b 

15 are shown in Figure 8b, but any number of contacts can be employed, with the provision 
that there are at least two contacts. The sticker can be connected to system 120 via a 
connector 143. Connector 143 is optionally foldable to facilitate packaging and storage 
of sticker 141. In various exemplary embodiments of the invention connector 143 
includes two conductive members 149a and 149b devoid of electrical communication 

20 therebetweeiL Each of electrical contacts 145a and 145b is in electrical communication 
with one conductive member of connector 143 via a difierent internal conducting line 
147a and 147b. Thus, in the present embodiments, sticker can be connected to system 
120 using a single line, because coimector 143 inter&ces communication for both 
contacts. 

25 Figure 8f is a schematic illustration of a package of several stickers (four such 

stickers are shown in the exemplary illustration of Figure 8f), where each sticker can be 
similar to sticker 141 described above. Also shown in Figure 8f are exemplary 
dimensions and distances of the package and the individual stickers. 

Following are technical preferred values which may be used for selective steps 
30 and parts of the embodiments described above. 

The output radiofrequency signal is preferably from about 10 KHz to about 
200 KHz in fiiequency and from about lOmV to about 200mV in magnitude; the input 
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radiofrequency signal is preferably about 75 KHz in frequency and about 20mV in 
magnitude; a typical impedance which can be measured by the present embodiments is 
from about 5 Ohms to about 75 Ohms; the resulting signal-to-noise ratio of flie present 
embodiments is at least 40dB; low pass filter 134 is preferably characterized by a cutoff 
frequency of about 35H2 and digitizer 138 preferably samples the signal at a rate of 
about 500-1000 samples per second. 

As used herem, the singular form "a", "an" and "the" include plural references 
unless the context clearly dictates otherwise. For example, the term "a compound" or "at 
least one compound" may include a plurality of compounds, includmg mixtures thereof 

Additional objects, advantages, and novel features of the present invention will 
become apparent to one ordinarily sldlled in the art upon examination of the foUowmg 
examples, which are not intended to be limiting. Additionally, each of the various 
embodiments and aspects of the present invention as delineated hereinabove and as 
cldmed in the claims section below finds experimental support in the follovrtng 
examples. 

EXAMPLES 

Reference is now made to the following examples, which together with the above 
descriptions, illustrate the invention m a non limiting fashion. 

EXAMPLE 1 
Prototype System 

A proto^e of a system for measuring blood flow in an organ of a subject 
according to the above description was constructed. 
The prototype system includes: 

(a) a self made radiofrequency generator generating an output 
radiofi^uency signal, TOKhz in frequency and 20mV in magnitude; 

(b) a plurality of electrodes, as described in Figures 8a-d; and 

(c) a double balanced mixer, purchased from Mini-Circuits, used for 
providing a radiofrequency sum and a radiofrequency difference, as detailed above. 

The prototype system fiirther includes electronic circuitry formed in a printed 
ckcuit board. Several electronic circuitries were designed and manufactured, so as to 
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investigate the correlation between the qualities of the results, the design of the 
electronic circuitry and the number of electrodes. The various electronic circuitries are 
schematically illustrated in Figures 9a-^. 

Figure 9a shows a block diagram of electronic circuitry to be used with three 
5 electrodes. The electrodes leads are designated in Figure 9a by E|, E2 and Ii, where the 
ou^ut radiofrequency signal, generated by the radiofirequency generator (designated 
OSC), is outputted through Ei and E2 and the input radiofrequency signal, as measured 
of the body is inputted through Ii. 

The input signal is channeled through a differential amplifier Gi, a band pass 

10 filter EPF and a differential amplifier G2. The input signal is channeled through a 
differential amplifier G3, a band pass filter EPF and an envelope elimination imit EEU. 
The EEU eliminates the amplitude modulation &om the input signal. Both input and 
output signals are mixed by mixer DMB, to form, as stated, a fiequency sum and a 
frequency difference. A low pass filter LPF filters out the frequency sum and the 

15 resulting signal (carrying the frequency difference) is further amplified by additional 
differential amplifiers G5, G^ and G7. Once amplified, tiie signal is digitized by an 
analog to digital digitizer and passed, via a USB conununication inter&ce to a 
processing and display unit. The processing unit includes a dynamically variable filter 
according to various exemplary embodiments of the present invention. 

20 Figure 9b shows a block diagram of electronic chrcuitry to be used with two 

electrodes of brain intra-luminal blood volume measurements. As there are only two 
electrodes E2 and Ii are combined to a single lead Ii. 

Hie ou^ut signal is channeled through a differential amplifier Gi, a band pass 
filter EPF and a differential amplifier G2. The input signal is channeled through a 

25 differential amplifier G2, a band pass filter EPF and an envelope elimination unit EEU 
which eliminates the amplitude modulation from the input signal. Both input and output 
signal are mixed by mixer DMB, to form the frequency sum and difference. The low 
pass filter LPF filters out the frequency sum and the resulting signal is fiirther amplified 
by additional differential amplifiers G4, G5 and Ge. As in the case of three electrodes, 

30 the signal is digitized by an analog to digital digitizer and passed, via a USB 
communication interface to a processing and display unit. 
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Figure 9c shows a block diagram of electronic circuitry to be used with four 
electrodes. The four leads, designated Ei, E2, Ii and I2, where the output radiofirequency 
signal, generated by radiofrequency generator OSC, is outputted through Ei and E2 and 
the input radiofiequency signal, as measured of the body are inputted through Ii and fc. 
In addition, the four leads, Ei, E2, Ii and I2 are connected to the body through capacitois 
designated C|, C2, C3 and C4. 

The principles of the cuxsuitry of Figure 9c are similar to the principles of the 
ckcuitry of Figure 9a with three electrodes. The advantage of the circuitry of Figure 9c 
is that by using both input leads Ii and I2 (as opposed to one input lead Ii of Figure 9a), 
effects of impedance differences between the electrodes and the body can be minimized 
Specifically, the influence of the voltage drop Ii and I2 is controlled by the characteristic 
impedance of the differential amplifier G3, which is selected to be sufGlciently large so 
that any impedance changes due to the contact between the body and the electrode is 
negligible, compared to the impedance of G3. 

Figure 9d shows a block diagram of the analog amplification circuit, which was 
used to amplify the radiofrequency signal after the low pass filtering in which the 
radiofrequency sum was filtered out 

EXAMPLE 2 
CUnical Trials 

The prototype system described in Example 1 was tested on human volunteers. 
The present Example includes a representative collection of trials performed on four of 
the volunteers. 
Methods 

Each subject was connected to four electrodes of the prototype system. Two 
electrodes served for input/output radiofrequency signals and two served as ECG leads. 

A radiofrequency signal pertaining to hemodynamic reactance was sampled at a 
sampling rate of 500 samples per second during continuous time intervals of 8 seconds. 
The signal was filtered by an analog low pass filter of 9 Hz. 

An ECG signal was sampled at the same rate (500 samples per seconds) and 
filtered using an analog filter of 250 Hz. 
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In all the trials, the signals acquired from each subject were filtered using two 
types of digital filters: a fixed filter with a lower bound of 0.9 Hz and an upper bound of 
6 Hz, and a dynamically variable filter in which the firequency bounds were varied in 
response to changes in the heart rate of the respective subject To this end, the linear 
5 dependence as illustrated in Figures 2a-b was used. 
Results 

Figures lOa-e show snapshots of the display of the prototype system obtained 
during a trial in which the electrodes of the system were connected to subject No. 1. 
Signals were acquired while the subject was stable (heart rate of 95 bpm). 

10 Figures lOa-b show results obtained using a fixed filter (Figure 10a) and 

dynaimcally variable filter (Figure 10b). In each of Figures lOa-b, there are seveii 
curves, designated, firom top to bottom, I, II, dl, dll, N, dN and ddN. The four top 
curves (1, 11, dl and dll) are ECO signals (leads I and II) and derivatives thereof (dl and 
dn, respectively). The three lowermost curves (N, dN and ddN) correspond to a 

15 hemodynamic reactance (N), its first time-derivative (dN) and second time-derivative 
(ddN). The right pane of Figure lOa show various calculated values, such as heart rate 
(BPM), cardiac output (CO), stroke volume (SV), ventricular ejection time (VETM), and 
the like. 

Figures lOc-d demonstrate the data analysis performed to provide the results 
20 presented in Figures lOa-b. Figure lOc demonstrates application of the fixed filter on 
the ECO signal. Figure lOd demonstrates the application of the fixed on the 
hemodynamic reactance signal, and Figure lOe demonstrates the application of the 
dynamically variable filter on the hemodynamic reactance signal. In each of Figures 
lOc-d there are five graphical representations, corresponding to, firom top to bottom: (i) 
25 the respective signal before filtering, (ii) the respective signal after filtering, (iii) the 
spectrum (Fourier decomposition) of the respective signal before filtering, (iv) the 
spectrum of the respective signal after filtering and (v) phase shift data. The lower and 
upper firequency bounds of the filter for the ECO signal were 1.2 Hz and 40 Hz, 
respectively; the lower and upper firequency bounds of the fixed filter for the 
30 hemodynamic reactance signal were 0.9 Hz and 6 Hz, respectively; and the lower and 
upper firequency bounds of the dynamically variable filter for the hemodynamic 
reactance signal were 1.4 Hz and 6.9 Hz, respectively. 
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In Figure 10c, the abscissae are scaled to 300 ms per division in the time domain 
repres^tations, and 3 Hz per division in all frequency domain representations. In 
Figures lOd-e, the abscissae are scaled to 300 ms per division in the time domain 
representations, and 0.5 Hz per division in the frequency domain representations. 
5 Figures lla-e show snapshots of the display of the prototype system obtained 

during another trial in which the electrodes of the ssystem were connected to subject No. 
1. In this trial, the subject was also stable (heart rate of 114 bpm). The graphical 
representations in Figures lla-e correspond to the same observables as Figures lOa-e. 

The lower and upper frequency bounds of the fixed filters were the same as in 
10 the trial above. The lower and upper frequency bounds of the dynamically variable filter 
for the hemodynamic reactance signal were 1.7 Hz and 7.4 Hz, respectively. 

Figures 12a-e show snapshots of the display of the prototype system obtained 
during another trial in which the electrodes of the system were connected to subject No. 
1. In this trial, the subject was agitated (heart rate of 140 bpm). The graphical 
15 representations in Figures 12a-e correspond to the same observables as Figuies lOa-e. 
The lower and upper frequency bounds of the fixed filters were the same as in the trials 
above. The lower and upper frequency bounds of the dynamically variable filter for the 
hemodynamic reactance signal were 2.1 Hz and 8 Hz, respectively. 

Figures 13a-e show snapshots of the display of the prototype system obtained 
20 during another trial in which the electrodes of the system were connected to subject No. 
1. In this trial, the subject was also agitated (heart rate of 137 bpm). The graphical 
representations in Figures 13a-e correspond to the same observables as Figures lOa-e. 
The lower and upper frequency bounds of the fixed filters were the same as in the trials 
above. The lower and upper fi^uency bounds of the dynamically variable filter for the 
25 hemodynamic reactance signal were 2.1 Hz and 7.9 Hz, respectively. 

Figures 14a-c show snapshots of the display of the prototype system obtained 
during a trial in which the electrodes of the system were connected to subject No. 2. In 
this trial, the subject was agitated (heart rate of 121 bpm). The graphical representations 
in Figures 14a-e correspond to the same observables as Figures lOa-e. The lower and 
30 upper frequency boimds of the fixed filters were the same as in the trials above. The 
lower and upper frequency bounds of the dynamically variable filter for the 
hemodynamic reactance signal were 1 .8 Hz and 7.5 Hz, respectively. 
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Figures 15a-e show snapshots of the display of the prototype system obtained 
during another trial in which die electrodes of the system were connected to subject No, 
2. In this trial, the subject was also agitated (heart rate of 123 bpm). The graphical 
representations in Figures 15a-e correspond to the same observables as Figures lOa-e. 
The lower and upper fiequenoy bounds of the fixed filters were the same as in the trials 
above. The lower and upper firequency bounds of the dynamically variable filter for the 
hemodynamic reactance signal were 1.8 Hz and 7.6 Hz, respectively. 

Figures 16a-g show snapshots of the display of the prototype system obtained 
during a trial in which the electrodes of the system were connected to subject No. 3. In 
this trial, the subject was agitated (heart rate of 121 bpm). The graphical repiesentations 
in Figures 16a-b correspond to the same observables as Figures lOa-b, and the graphical 
representations in Figures 16e-g correspond to the same observables as Figures lOc-e. 
Figures 16c-d are respective zoom-in images of Figures 16a-b. The lower and upper 
frequency bounds of the fixed filters were the same as in the trials above. The lower and 
upper firequency bounds of the dynamically variable filter for the hemodynamic 
reactance signal were 1.8 Hz and 7.5 Hz, respectively. 

Figures 17a-g show snapshots of the display of the prototype system obtained 
during a trial in which the electrodes of the system were connected to subject No. 4. In 
this trial, the subject was agitated (heart rate of 139 bpm). The graphical representations 
in Figures 17a-g correspond to the same observables as Figures 16a-g. Figures 17c-d are 
respective zoom-in images of Figures 17-b- The lower and upper frequency bounds of 
the fixed filters were the same as in the trials above. The lower and upper fi^uency 
bounds of the dynamically variable filter for the hemodynamic reactance signal were 
2.1 Hz and 8 Hz, respectively. 

As shown in Figures I0a-17g, the dynamically varying filter significantly 
improves the qiiality of the results, particularly when the subjects are agitated (Figures 
12a- 17g). Additionally, the dynamically variable filtering technique of the present 
embodiments allows consistent calculation of CO values. For example, as demonstrated 
in the trials with Subject No. 3 (see Figures 16a-e), the two filtering techniques resulted 
in different CO values: 35.6 L/min for the fixed fdtering technique and 21.99 L/min for 
dynamically variable filtering technique. Similar improvements were observed in other 
subjects. 
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Plots of cardiac output as calculated from signals filtered using the two filtering 
techniques are shown in Figures 18a-b, 19a-b and 20a-b for subjects Nos. 2, 3 and 4, 
respectively. Figures 18a, 19a and 20a show results obtained using the dynanucally 
variable filter of the present embodiments and Figures 18b, 19b and 20b show results 
obtained using fixed filter. As demonstrated, the CO values obtained using the 
dynamically variable filter of the present embodiments are more accurate. 

Table 1 below summarizes the cardiac output results obtained for the above 

trials. 



Table 1 



Trial No. 


Subject 


Heart Rate 
[bpm] 


filter 


CO [L/min] 


type 


fiequency band \JHz] 


1 




95 


fixed 


0.6-9.0 


14.62 


variable 


1.4-6.9 


14.68 


2 




114 


fixed 


0.6-9.0 


18.81 


variable 


1.7-7.4 


19.77 


3 




140 


fixed 


0.6-9.0 


18.55 


variable 


2.1-8 


18.41 


4 




137 


fixed 


0.6-9.0 


15.82 


variable 


2.1-7.9 


17.05 


5 


2 


121 


fixed 


0.6-9.0 


13.92 


variable 


1.8-7.5 


18.47 


6 


2 


123 


fixed 


0.6-9.0 


19.42 


variable 


1.8-7.6 


17.84 


7 


3 


121 


fixed 


0.6-9.0 


35.60 


variable 


1.8-7.5 


21.99 


8 


4 


139 


fixed 


0.6-9.0 


24.94 


variable 


2.1-8.0 


24.37 



Although the invention has been described in conjunction with specific 
embodiments thereof, it is evident that many altematives, modifications and variations 
will be apparent to those skilled m the art Accordingly, it is intended to embrace all 
such altematives, modifications and variations that fell within the spirit and broad scope 
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of the appended claims. All publications, patents and patent applications mentioned in 
this specification are herein incorporated m their entirety by reference into the 
specification, to the same extent as if each individual publication, patent or patent 
application was specifically and individually indicated to be incorporated herein by 
5 reference. In addition^ citation or identification of any reference in this application shall 
not be construed as an admission that such reference is available as prior art to the 
present invention. 
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WHAT IS CLAIMED IS: 

1. A method of processing an input signal pertaining ^to at least one 
electrical property of an organ of a subject, comprising determining a physiological 
condition of the subject, selecting a frequency band, filtering said signal according to 
said frequency band, and dynamically adapting said frequency band in response to a 
change in said physiological condition, thereby processing the signal. 

2. A filtering device, comprising: a first input unit for receiving an input 
pertaining to at least one electrical property of an organ of a subject, a second input unit 
for receiving data pertaining to a physiological concfition of the subject, and a filtering 
unit configured for filtering said input signal according to a frequency band which is 
djnnamically adapted in response to a change in said physiological condition. 

3 . A system for monitoring cardiac output, comprising the device of claim 2. 

4. A system for predicting at least one of: a body cell mass, a frit free mass 
and total body water of a subject, comprising the device of claim 2. 

5. A system for determining hematocrit of blood in a body part of a subject, 
comprising the device of claim 2. 

6. A system for monitoring hydration status of a subject, comprising the 
device of claim 2. 

7. A system for discriminating tissue, comprising the device of claim 2. 

8. A system for calculating the circimiference of a body segment, 
comprising the device of claim 2. 

9. A method of monitoring at least one electrical property of an organ of a 
subject, comprising sensing an input radiofrequency signal from the organ, processing 
said input radiofrequency signal to provide a processed input signal, filtering said input 
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signal using a dynamically variable filter to provide a filtered signal, and using said 
filtered signal for monitoring the at least one electrical property of the organ. 

10. The method of claim 9, wherein said dynamically variable filter is 
adapted in response to a change in a physiological condition of the subject 

1 1. The method of claim 9> wherein said possessing comprises filtering said 
input radiofirequency signal using an analog filter. 

12. The method of claim 9, fiirther comprising calculating at least one 
quantity^ using said electrical property, wherein said at least one quantity is selected from 
the group consisting of a stroke volume, a cardiac output, a brain intra luminal blood 
volume and blood flow. 

13. The method of clmm 9, wherein said blood flow comprises at least one 
of: an external carotid blood flow rate, an internal carotid blood flow rate, an ulnar blood 
flow rate, a radial blood flow rate, a brachial blood flow rate, a common iliac blood flow 
rate, an external iliac blood flow rate, a posterior tibial blood flow rate, an anterior tibial 
blood flow rate, a peroneal blood flow rate, a lateral plantar blood flow rate, a medial 
plantar blood flow rate and a deep plantar blood flow rate. 

14. The method of claim 12, fiirther comprising detennining a phase shift of 
said input radiofirequency signal relative to an output radiofirequency signal transmitted 
to the organ and using said phase shift for calculating said at least one quantity. 

15. The method of claim 14, fiirther con^)rising transmitting said output 
radiofrequency signal to the organ. 

16. The method of claim 14, fiulher comprising reducing or eliminating 
amplitude modulation of said input radiofirequency signal, so as to provide a signal of 
substantially constant envelope. 
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17. The method of claim 16, wherein said reducing or eliminating said 
amplitude modulation comprises maintaining a phase modulation of said input 
radiofrequency signal. 

1 8. The method of claim 9, further comprising nuxing said input 
radiofrequency signal and an output mdiofrequency signal transmitted to the organ so as 
to provide a mixed radiofrequency signaL 

19. The method of claim 18, wherein said mixing comprises providing a 
radiofrequency sum and a radiofrequency difference. 

20. Apparatus for monitoring at least one electripal property of an organ of a 
subject, comprising: 

an input unit for receiving an input radiofrequency signal sensed from the organ; 

a signal processing unit for processing said input radiofrequency signal to 
provide a processed input signal; 

a filtering unit configured for filtering said input signal using a dynamically 
variable filter to thereby provide a filtered signal; and 

a monitoring unit for monitoring the at least one electrical property of the organ 
based on said filtered signal. 

21- The apparatus of claim 20, wherein smd dynamically variable filter is 
adapted in response to a change in a physiological condition of the subject 

22. The method, device, system or apparatus of any of claims 1, 2-8, 10 and 
21, wherein said physiological condition is a heart rate of the subject. 

23. The method, device, system or apparatus of claim 22, wherein at least one 
of a lower frequency bound and an upper frequency bound characterizing said filter is a 
linear frinction of said heart rate. 
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24. The method, device, system or apparatus of claim 23, wherein said lower 
frequency bound is about 0-9*(HR/60) Hz at all times, wherein said HR is said heart rate 
in units of beats per minute. 

25. The method, device, system or apparatus of claim 23, wherein s^d upper 
frequency bound is about 6 + l.S*[(HR/60) - 1] Hz at all times, wherein said HR is said 
heart rate in units of beats per minute. 

26. The method, device, system or apparatus of claim 22, wherein said heart 
rate is determined from an ECO signal received from said subject. 

27. The method, device, system or apparatus of any of claims 1, 2-8, 10 and 
21, wherein an upper frequency bound characterizing said filter is determined using an 
iterative process. 

28. The method, device, system or apparatus of claim 27, wherein said 
iterative process is based on a comparison between a value of a physiological parameter 
as extracted from said filtered input signal and a value of said physiological parameter as 
extracted from a reference signal. 

29. The method, device, system or apparatus of claim 28, wherein said 
reference signal comprises an ECG signal. 

30. The method, device, system or apparatus of claim 27, wherein each 
itemtion of said iterative process comprises: 

if said comparison meets a predetermined criterion , then updating said upper 
frequency bound by calculating an average between a low threshold and a high 
threshold; 

said thresholds being predetermined or set in a previous iteration of said iterative 
process. 
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31. The method, device, system or apparatus of claim 28, wherein said 
physiological parameter is a ventricular ejection time (VET). 

32. The method, device, system or apparatus of claim 28, wherein said VET 
is averaged over a plurality of heart beats. 

33. The method, device, system or apparatus of claim 28, wherein said VET 
is extracted from an average heart beat morphology of the subject 

34. The method, device, system or apparatus of claim 27, wdierein said 
physiological condition is a heart rate of the subject 

35. The method, device, system or apparatus of claim 34, wherein an initial 
value of said upper frequency bound in said iterative process is a linear function of said 
heart rate. 

36. The method, device, system or apparatus of claim 35, wherein an initial 
value of said upper frequency bound in said iterative process is about 6 + L5*[(HR/60) - 
1] Hz, wherein said HR is said heart rate in units of beats per minute. 

37. The mefliod or apparatus of clahn 9 or 20, v^ierem said input 
radiofiequency signal is indicative of impedance the organ. 

38. The method or apparatus of claim 9 or 20, vsdierein said input 
radiofiequency signal is indicative of hemodynamic reactance of the organ. 

39. The apparatus of claim 20, wherein said signal possessing unit comprises 
an analog filter for filtering said input radiofrequency signal. 

40. The apparatus of claim 20, further comprising a data processor for 
calculating at least one quantity usmg said electrical property, wherein said at least one 
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quantity is selected from the group consisting of a stroke volume, a cardiac ou^ut, a 
brain intra limiinal blood volume and blood flow. 

41. The apparatus of claim 40, further compri^g a phase shift determinator 
configured for determining a phase shift of said input radiofrequency signal relative to 
an output radiofrequency signal transmitted to the organ, herein said data processor is 
configured for calculating said at least one quantity based on said phase shift 

42. The apparatus of claim 20, wherein said signal processing imit comprises 
an envelope elimination unit designed and configured to reduce or eliminate amplitude 
modulation of said filtered signal, so as to provide a signal of substantially constant 
envelope. 

43. The apparatus of claim 42, wherein said envelope elimination unit is 
designed and configured to maintain a phase modulation of said signal. 

44. The appamtus of claim 20, wherein said signal processing unit comprises 
a mixer designed and configured to mix said input radiofrequency signal and said output 
radiofrequency signal transmitted to the organ so as to provide a mixed radiofrequency 
signal. 

45. The apparatus of claim 44, wherein said mixer is operable to provide a 
radiofrequency sum and a radiofrequency difference. 

46. A system for monitoring at least one electrical property of an organ of a 
subject, comprising: 

a radiofrequency generator for generating an ouQ)ut radiofrequency signal; 

a plurality of electrodes, designed to be comiectable to the skin of the subject, 
said electrodes being for transmitting said output radiofrequency signal to the organ and 
for sensing an input radiofrequency signal from the organ; and 

the apparatus according to any of claims 20-45. 
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